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ABSTRACT 


This  report  covers  in  detail  the  solid  state  research  work 
at  Lincoln  Laboratory  for  the  period  1  October  1964  through 
31  December  1964.  The  topics  covered  are  Solid  State  De¬ 
vice  Research,  Laser  Research,  Materials  Research,  Band 
Structure  and  Spectroscopy  of  Solids,  and  Magnetism  and 
Resonance. 
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INTRODUCTION 


I.  SOLID  STATE  DEVICE  RESEARCH 

Coherent  emission  at  5.2  microns  parallel  to  the  direction  of  the  diode  current  has  been 
obtained  in  n+pp+  InSb  diodes.  Large-volume  luminescence  and  coherent  emission  in  similar 
n+pp+  InSb  diodes  were  described  in  the  last  quarterly  progress  report,  but  the  coherent  beam 
was  emitted  perpendicular  to  the  direction  of  current,  as  in  all  other  existing  semiconductor 
lasers.  However,  the  large  active  regions  of  these  InSb  diodes,  in  which  lifetimes  of  10  ^  to 
10  6  sec  permit  population  inversion  in  regions  of  the  order  of  100  microns  away  from  the  in¬ 
jecting  contacts,  are  naturally  suited  for  laser  oscillations  longitudinal  to  the  current.  Such 
lasers  are  free  from  some  of  the  structural  restrictions  of  the  perpendicular  lasers;  for  ex¬ 
ample,  they  can  be  more  readily  constructed  in  arrays.  They  also  promise  to  emit  coher¬ 
ently  over  large  areas,  hence  with  small  beam  angles,  and  they  should  be  suited  for  large 

power  outputs.  Coherent  emission  was  obtained  near  10°K  with  50-nsec  current  pulses  of 
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20  amp  (6  x  10  amp  cm  ),  with  a  minimum  longitudinal  magnetic  field  of  7  kgauss.  Com¬ 
parison  of  the  wavelength  and  the  magnetic  shift  with  earlier  data  on  InSb  diodes  indicates  that 
the  highest  energy  state  of  the  spin- split  lowest  Landau  level  is  involved  in  the  transition. 

o 

Mode  structure  was  observed  with  mode  spacing  of  120  A  as  expected  from  the  220-micron 
length  of  the  longitudinal  cavity. 

The  PbSe  diode  laser,  with  a  coherent  emission  peak  at  8.5  microns  at  12°K,  has  been  further 

developed.  Fabrication  techniques  are  based  on  controlling  carrier  type  and  concentration  by 

adjusting  the  Pb:Se  ratio,  rather  than  adding  impurities.  The  diode  emission  below  threshold 

exhibits  two  spectral  peaks,  one  with  a  maximum  near  8.5  microns  which  increases  superlin- 

early  with  current,  and  another  with  a  maximum  at  10. 1  microns  which  increases  slowly  with 

current.  Laser  action  associated  with  the  8.5-micron  peak  occurs  above  a  threshold  current 
_2 

density  of  2000  amp  cm  .  The  unresolved  emission  peak  shifts  linearly  to  higher  energies 

-8 

with  a  magnetic  field  oriented  in  the  <100>  direction  at  the  rate  of  7. 1  x  10  ev/gauss,  or 
17  Mcps/gauss,  which  is  the  expected  value  if  the  emission  is  associated  with  band-to-band 
transitions.  The  threshold  current  density  decreases  to  a  fraction  of  its  zero  field  value  in  a 
magnetic  field  of  approximately  10  kgauss,  then  rises  slowly  with  increasing  field.  The  PbSe 
laser  should  eventually  be  of  particular  interest  for  communications  since  its  emission  is  in 
the  8-  to  14-micron  atmospheric  window,  a  spectral  region  of  high  atmospheric  transparency 
where  attenuation  due  to  scattering  by  haze  is  generally  low. 


II.  LASER  RESEARCH 

The  frequency,  far-field,  and  intensity  characteristics  of  stimulated  Brillouin  emission  at 
90®  to  the  ruby  laser  beam  in  nitrobenzene  have  been  measured.  The  observed  Brillouin  fre¬ 
quency  shift  is  consistent  with  that  measured  for  90°  Brillouin  scattering  of  a  He-Ne  laser  in 
nitrobenzene. 


Further  etalon  studies  of  the  stimulated  Raman  emission  at  90°  have  shown  the  presence  of 
many  on-axis  Raman  cavity  modes  as  well  as  a  continuum  of  walkoff  modes,  all  within  a  few 
cm  bandwidth. 

Theoretical  studies  of  the  stimulated  Stokes  and  anti-Stokes  radiation  from  a  point  source  pre¬ 
dict  exponential  growth  of  a  spherical  Stokes  wave,  and  exponential  damping  of  anti-Stokes  ra¬ 
diation  in  all  directions  including  that  of  phase  matching.  However,  the  anti- Stokes  radiation 
will  be  emitted  near  a  boundary  of  the  Raman  medium. 

Two-photon  absorption  effects  were  observed  in  the  photoelectron  current  of  a  multiplier  illu¬ 
minated  by  a  focused  CW  infrared  He-Ne  laser.  Unusual  photoconduction  effects,  partly  at¬ 
tributable  to  two-photon  absorption,  were  observed  in  several  samples  of  large-gap,  high- 
resistivity  materials  illuminated  by  a  spiking  ruby  laser. 

Measurements  of  the  amplitude  noise  of  a  stable  He-Ne  laser  have  shown  good  agreement  with 
bandwidth  and  noise  intensity  predictions  of  a  van  der  Pol  oscillator  theory  for  above  threshold 
data  and  a  linear  maser  amplifier  theory  for  below  threshold  data. 

Stable  pulsed  arc  discharge  conditions  have  been  achieved  in  the  high-pressure  mercury  sys¬ 
tem  to  obtain  an  estimated  single  pass  gain  of  0.  3  percent.  A  sufficiently  low  loss  optical 
cavity  is  being  assembled  for  obtaining  laser  oscillation. 

III.  MATERIALS  RESEARCH 

Single  crystals  of  sapphire  (melting  point  2050° C)  have  been  grown  by  the  Czochralski  method 
in  a  resistance  furnace  that  was  recently  developed  for  operation  in  oxidizing,  neutral,  or  re¬ 
ducing  atmospheres  at  temperatures  up  to  2400° C.  The  melt  was  contained  in  an  iridium  or 
molybdenum  crucible  under  an  argon  atmosphere.  Since  the  molten  sapphire  did  not  evaporate 
appreciably,  it  was  possible  to  remove  the  zirconia  tube  used  to  protect  the  tantalum  heater 
element  from  oxygen  and  other  reactive  gases. 

Additional  experiments  on  the  vapor  growth  of  iodine  crystals  by  forced  convection  have  pro¬ 
vided  further  confirmation  that  constitutional  supercooling  of  the  vapor  occurs  at  sufficiently 
high  growth  velocities.  Additional  data  have  been  obtained  on  the  variation  of  growth  velocity 
with  the  temperature  difference  between  the  source  of  iodine  vapor  and  the  growing  interface. 
These  data  are  fitted  satisfactorily  by  an  equation  derived  from  elementary  boundary  layer 
theory. 

3  3 

Single  crystals  of  T^O^  (~1  cm  )  and  TiO  (~1  mm  )  have  been  grown  by  the  Czochralski 
method  in  the  resistance  furnace  used  for  sapphire  growth.  Large-grained  boules  of  Ti^O^ 
have  been  obtained  by  the  Vemeuil  method  with  an  induction  plasma  torch.  Measurements  of 
the  Seebeck  coefficient  have  been  made  on  a  single  crystal  of  T^O^  between  7°  and  550° K. 

Hall  coefficient  (R^)  measurements  have  been  made  at  77 °K  on  several  samples  of  In^SbTe^. 

Application  of  the  one-carrier  expression  n  =  -1/R„e  gives  electron  concentrations  of  1. 1  and 

22  -3  ^  21  -3 

1.2  x  10  cm  for  two  high-pressure  samples  and  8.8,  9.0  and  9.  4  x  10  cm  for  three 

low-pressure  samples.  These  results  indicate  that  the  difference  in  superconducting  behavior 

between  high-  and  low-pressure  samples  may  be  due  to  a  difference  in  carrier  concentrations. 


vi 


The  partial  pressures  of  Hg(g)  and  Te2(g)  in  equilibrium  with  Hg-saturated  HgTe(c),  Te- 
saturated  HgTe(c),  and  someHg-Te  melts  between  about  400°  and  700°C  have  been  determined 
by  measuring  the  optical  density  of  the  vapor  as  a  function  of  wavelength.  The  partial  pres¬ 
sures  depend  strongly  on  the  composition  of  the  solid  phase.  The  Gibbs  free  energy  of  forma¬ 
tion  of  HgTe(c)  from  Hg(g)  and  Te2(g)  is  found  to  be  AG°  =-41.66  +  42. 71(10” ^)Tkcal/mode. 

Wet  chemical  methods  have  been  developed  for  determining  Hg  and  Te  in  HgTe  with  an  accu¬ 
racy  of  about  1  part  per  thousand  and  for  determining  In,  Sb  and  Te  in  0. 2  g  In-Sb-Te  samples 
with  an  accuracy  of  about  2  parts  per  thousand  for  each  component. 


IV.  BAND  STRUCTURE  AND  SPECTROSCOPY  OF  SOLIDS 

The  study  of  the  interband  transitions  in  antimony  initially  reported  in  the  last  quarter  has 
been  further  extended.  The  most  probable  locations  in  momentum  space  of  the  energy  minima 
involved  in  the  interband  transitions  have  been  deduced  and  the  reduced  effective  masses  and 
energy  gaps  have  been  obtained.  It  has  also  been  found  that  the  energy  bands  are  nonparabolic. 
In  order  to  more  fully  appreciate  these  results,  the  dispersion  relations  for  the  relevant  points 
in  the  Brillouin  zone  have  been  developed  in  terms  of  the  k  •  p  Brillouin-Wigner  perturbation 
theory  and  the  symmetry  properties.  The  nonparabolic  character  of  the  energy  bands  is  in¬ 
corporated  in  this  development;  furthermore,  it  is  possible  to  include  the  effects  of  spin-orbit 
interaction  in  this  formalism. 

In  another  related  problem,  the  magnetoplasma  effects  in  bismuth  are  being  investigated  in 
detail.  Studies  on  this  semimetal  offer  the  advantage  that  the  contributions  due  to  interband 
effects  should  be  negligible  near  the  plasma  frequencies.  Thus,  this  work  should  serve  as  a 
guide  in  interpreting  the  more  complex  effects  observed  in  antimony  where  the  interband  and 
intraband  effects  are  coupled  to  each  other. 

The  Faraday  rotation  in  the  energy  region  of  the  indirect  transition  in  germanium  is  being 
studied  at  high  magnetic  fields;  oscillatory  effects  due  to  the  exciton  and  also  due  to  transi¬ 
tions  from  the  valence  band  to  the  electronic  Landau  ladder  of  the  conduction  band  have  been 
observed.  Since  Faraday  rotation  is  a  dispersive  effect,  contributions  due  to  the  direct  and 
indirect  transitions  are  superimposed  and  thus  a  direct  comparison  of  the  matrix  elements  in¬ 
volved  in  the  respective  transitions  can  be  made. 

In  GaSb  the  study  of  the  magnetoabsorption  due  to  exciton  complexes,  believed  to  consist  of  an 
ionized  acceptor  ion  and  an  exciton,  has  been  extended  to  90  kgauss.  It  has  been  found  that  the 
functional  behavior  of  the  diamagnetic  shift  changes  at  magnetic  fields  above  ~20  kgauss  and 
tends  to  a  linear  behavior;  this  is  in  agreement  with  the  theory  of  Elliott  and  Loudon.  In  a  re¬ 
lated  theoretical  development  of  the  optical  absorption  due  to  excitons,  the  theory  of  Elliott, 
which  holds  in  the  limit  of  small  wave  vector,  has  been  extended  to  the  more  general  case  of 
finite  wave  vector. 

The  magnetic  field  dependence  of  the  photoluminescence  of  GaAs  has  been  studied  for  both  n- 
and  p-type  material.  In  n-type  material  the  exciton  lines,  as  well  as  the  transitions  proceed¬ 
ing  through  holes  trapped  at  acceptors,  exhibited  a  quadratic  behavior  at  low  fields  which 
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eventually  became  linear  at  high  fields,  and  intensities  which  were  independent  of  magnetic 
field.  This  type  of  energy  shift  can  be  accounted  for  by  a  simple  effective  mass  hydrogenic 
model;  however,  hole  spin  splittings  which  were  expected  on  the  basis  of  previous  experiments 
were  not  observed.  In  p-type  materials  an  additional  strong  line  was  observed;  as  the  mag¬ 
netic  field  was  increased  this  line  shrank  in  intensity,  while  the  other  lines  grew,  sometimes 
by  as  much  as  one  order  of  magnitude. 

In  an  attempt  to  resolve  the  polaron  anomaly  in  CdS,  a  Frohlich-type  Hamiltonian  with  an  iso¬ 
tropic  electron-phonon  interaction  term  is  being  used  as  a  model  to  calculate  the  piezoelectric 
electron-phonon  coupling.  The  effects  of  such  coupling  on  the  electron  energy  are  being  in¬ 
vestigated  at  zero  temperature  where  the  lowest  order  perturbation  theory  can  be  rigorously 
evaluated  unlike  the  Hopfield-Mahan  calculation  which  required  severe  approximations.  Very 
preliminary  results  indicate  that,  at  magnetic  field  strengths  of  several  kgauss,  this  model 
yields  a  piezoelectric  interaction  which  is  much  too  weak  to  account  for  the  experimental 
results. 

In  connection  with  the  band  structure  of  the  rare  earth  metals,  a  formal  solution  of  the  rela¬ 
tivistic  Schrddinger  equation  has  been  obtained  for  a  periodic  "muffin  tin"  potential,  using  the 
augmented  plane-wave  technique.  The  problem  does  not  appear  to  be  much  more  difficult  than 
that  of  evaluating  the  nonrelativistic  problem.  The  principal  source  of  inaccuracies  remains 
the  same,  namely,  the  lack  of  precision  in  the  assumed  potential.  Since  this  technique  ap¬ 
pears  feasible,  the  programming  of  the  calculation  to  obtain  the  energy  bands  of  specific  ma¬ 
terials  is  now  being  considered. 

The  study  of  deep  impurities,  such  as  sulfur,  in  silicon  is  continuing.  The  behavior  under 
stress  of  the  transitions  between  the  ground  state  of  an  impurity  center  and  its  excited  p-like 
states  has  been  examined.  It  has  been  found  that  all  transitions  behave  identically  under  stress 
and  that  this  is  consistent  with  a  simple  deformation  potential  model.  An  analysis  of  the  re¬ 
sults  has  yielded  a  pure  shear  deformation  potential  of  7.  9  ±  0.  2  ev. 

Galvanomagnetic  measurements  have  been  made  on  newly  grown  HgTe.  In  some  samples,  at 
4°K,  Hall  mobilities  larger  than  those  found  in  InSb  have  been  measured.  Both  n-  and  p-type 
material  have  been  studied;  magnetoresistance,  magneto-Hall  effects,  and  the  electron  Hall 
mobility  have  been  measured  at  4°  and  78°  K  in  fields  up  to  105  kgauss.  Magneto  re  si  stance 
ratios  of  up  to  200  have  been  observed.  Attempts  are  being  made  to  fit  the  data  to  a  mixed 
conduction  model.  The  Hall  coefficient  as  a  function  of  temperature  has  been  employed  to  ob¬ 
tain  a  value  for  the  thermal  energy  gap  Ef  in  HgTe,  using  an  analysis  similar  to  that  used  in 
HgSe.  The  value  obtained  is  ---0.02  ±  0.01  ev. 

V.  MAGNETISM  AND  RESONANCE 

Calculations  of  the  Curie  temperature  in  cubic  spinels  by  means  of  a  new  high-temperature 
expansion  have  been  completed  and  compared  with  results  by  older  methods.  The  improvement 
is  marked  for  the  case  of  zeroB-B  interactions,  but  is  pronounced  for  nonzero  B-B  interactions. 
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It  has  been  pointed  out  that  the  Overhauser- Stearns  conclusion  from  their  Mossbauer  meas¬ 
urements  on  Fe-Al  alloys  is  premature  because  the  plane-wave  approximation  for  the  conduc¬ 
tion  electrons  is  not  good  in  the  vicinity  of  the  nucleus.  If  the  plane  wave  is  orthogonalized  to 
the  core  electrons  for  determining  the  value  of  the  wave  function  at  the  nuclear  positions,  the 
spin  density  of  the  conduction  electrons  that  is  induced  by  the  localized  electrons  is  in  much 
better  accord  with  experiment.  It  was  shown  qualitatively  that  corrections  to  the  other  approx¬ 
imations  in  the  Ruderman-Kittel  theory  also  improve  the  agreement. 

The  need  for  improved  materials  for  digital  phase  shifters  for  radar-antenna  arrays  has  moti¬ 
vated  a  reinvestigation  of  the  conditions  for  achieving  spontaneously  square  B-H  loops  in  fer¬ 
rites.  It  was  shown  that  reverse-domain  nucleation  may  occur  either  at  a  grain  boundary  or 
at  a  chemical  inhomogeneity  having  roughly  0. 1  a  grain  volume,  a  magnetization  about  2  per¬ 
cent  different  from  the  parent  matrix,  and  no  static  crystallographic  discontinuity  at  the  inho¬ 
mogeneity  matrix  interface.  Evidence  has  been  found  for  the  existence  of  such  chemical  in¬ 
homogeneities  in  square-loop  ferrites  containing  small  concentrations  of  Jahn-Teller  ions. 

Susceptibility,  anisotropy,  and  resonance  data  for  the  ilmenites  FeTiO^  and  CoTiO^  have  been 
reconciled  by  the  assumption  that  the  spin-orbit  splitting  is  larger  than  the  trigonal -field 
splitting. 

Single  crystals  of  ReOQ  have  been  prepared  and  found  to  be  metallic  with  a  room-temperature 
—6  ^ 

resistivity  of  8  x  10  u  ohm-cm. 

The  perovskite  LaCoO^  was  found  to  exhibit  a  first-order  rhombohedral  ^  rhombohedral  phase 
change  at  1210°K,  and  the  crystallographic  data  appear  to  support  the  creation  of  two  cobalt 
sublattices  below  1210°K,  on  one  of  which  the  cobalt  ions  are  primarily  in  a  low-spin  state  and 
on  the  other  in  a  high-spin  state.  The  details  of  the  low-spin  ^  high-spin  distribution  as  a 
function  of  temperature  are  not  yet  firmly  established. 

Crystallographic  and  magnetic  studies  of  the  systems  FexMn^_xP  and  FeP-FeAs-FeSb  have 
been  initiated.  In  the  first  system  the  metamagnetic  ^  ferromagnetic  transition  increases 
with  x  to  reach  the  Curie  temperature  at  about  x  =  0. 25.  Within  the  range  0.  8  <  x  <  1.0  there 
is  a  change  from  metamagnetism  to  ferromagnetism. 

Two  main  areas  are  being  investigated  for  the  improvement  of  measurements  using  phonon 
generation  at  70  Gcps:  thin-film  transducers  and  the  reduction  of  receiver  noise.  The  system 
is  also  being  modified  to  permit  measurements  down  to  1.5°K. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  LONGITUDINAL  INJECTION-PLASMA  LASER  OF  InSb 

We  have  obtained  coherent  emission  in  n+pp+  InSb  structures  parallel  to  the  direction  of  the 
diode  current,  as  illustrated  in  Fig.I-1.  Large-volume  luminescence  and  coherent  emission 
were  previously  obtained  in  similar  n+pp+  InSb  structures,  but  the  coherent  beam  was  emitted 

i 

perpendicular  to  the  direction  of  current,  as  in  all  other  existing  semiconductor  lasers.  How- 
ever,  the  large  active  regions  of  these  InSb  diodes  in  which  lifetimes  of  10  to  10 ~  sec  permit 
population  inversion  in  regions  of  the  order  of  100  microns  away  from  the  injecting  contacts  are 
naturally  suited  for  laser  oscillations  longitudinal  to  the  current,  provided  a  suitable  Fabry- 
Perot  cavity  can  be  constructed  in  that  direction  and  provided  lasing  in  the  perpendicular  direc¬ 
tion  is  inhibited.  Such  lasers  are  free  from  some  of  the  structural  restrictions  of  the  perpen¬ 
dicular  lasers.  For  example,  they  can  be  more  readily  constructed  in  arrays.  They  also 
promise  to  emit  coherently  over  large  areas,  and  with  small  beam  angles,  and  they  should  be 
suited  for  large  power  outputs. 


Fig.  1-1.  InSb  n  pp  structure  with 
Fabry- Perot  cavity  designed  for  co¬ 
herent  emission  parallel  to  current. 


As  shown  in  Fig.  1-1,  the  laser  is  mounted  on  a  copper  heat  sink  over  a  small  hole.  An  n+ 

13  -3 

layer  is  grown  on  a  (100)  plane  of  the  p-type  base  (N^  —  N^  ~  3  X  10  cm  at  77°K)  by  a  solution- 
regrowth  technique  similar  to  that  developed  by  Nelson^  for  GaAs  and  germanium.  A  p+  layer 
is  formed  by  zinc  diffusion.  The  laser  cavity  consists  of  an  optically  flat  polished  surface  on 
the  n+  side  and  a  parallel  surface  on  the  p+  side  polished  prior  to  the  zinc  diffusion.  The  pf 
surface  is  made  totally  reflecting  by  evaporating  layers  of  silver  and  gold  on  a  small  circle, 
which  also  serves  as  an  electrical  contact.  Contact  to  the  n+  layer  is  made  by  soldering  directly 
to  the  copper  heat  sink  with  an  In-Te  alloy.  To  reduce  the  transverse  dimension  of  the  active 
region,  the  p+  layer  outside  the  evaporated  circle  is  removed  by  etching.  In  addition,  the  follow¬ 
ing  factors  tend  to  prevent  radial  spreading  of  the  plasma,  thereby  producing  an  active  region 
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which  is  long  in  the  direction  of  the  current,  extending  between  the  p+  circle  and  the  high  con- 
18  -3  -f 

ductivity  (Nn  «  7  X  10  cm  )  n  region.  First,  the  lifetime  of  electrons  inside  the  plasma  is 
3  U 

about  10  times  shorter  than  outside,  which  tends  to  localize  the  injected  carriers  as  noted 
earlier  in  such  plasmas.  A  second  collimating  effect  is  due  to  the  external  magnetic  field  which 
was  applied  longitudinal  to  the  plasma,  as  indicated  in  Fig.I-1.  The  presence  of  this  magnetic 
focusing  was  observed  in  the  voltage-current  characteristics  of  the  devices  with  fields  of  a  few 

4-6 

kilogauss.  A  third  factor  is  the  pinch-effect  of  the  azimuthal  self-magnetic  field  in  the  plasma. 
Pinching  of  this  kind,  however,  is  counteracted  by  the  applied  longitudinal  magnetic  field,  and 
does  not  appear  to  be  a  major  factor  under  conditions  existing  in  these  devices,  particularly 
since  appreciable  pinching  would  prevent  lasing  in  the  "transverse"  lasers. 


Fig.  1-2.  Mode  structure  of  longitudinal  laser 
of  Fig.  1-1  for  three  values  of  magnetic  field 
applied  parallel  to  current.  Laser  was  operated 
near  10°K  with  45-amp,  50-nsec  pulses. 


Coherent  emission  was  obtained  near  10°K  with  50-nsec  current  pulses  of  20  amp  (6  X 
4  -2 

10  amp  cm  ),  with  a  minimum  longitudinal  magnetic  field  of  7kgauss.  Figure  1-2  shows  the 

mode  structure  for  three  values  of  magnetic  field.  Comparison  of  the  wavelength  and  the  mag- 

7  8 

netic  shift  with  earlier  data  on  InSb  diodes  '  indicates  that  the  highest  energy  state  of  the  spin- 
split  lowest  Landau  level  is  involved  in  the  transition.  The  mode  spacing  near  13kgauss  is 

o  o 

120  A  and  remains  very  nearly  at  this  value  as  the  wavelength  is  shifted  over  1300  A  by  a  magnetic 

9 

field  change  of  11.3kgauss.  As  in  earlier  experiments,  the  mode  position  changed  slightly  with 
magnetic  field  due  to  a  change  in  the  index  of  refraction.  For  the  longitudinal  cavity  length  of 
220  microns,  the  experimental  mode  spacing  corresponds  to  a  value  of  5.2  for  [nQ  —  \Q(dn/dX)]. 
The  good  agreement  of  this  value  with  values  for  [nQ  —  XQ(dn/d\)]  obtained  in  transverse  InSb 
lasers  allows  us  to  conclude  that  oscillations  indeed  occur  in  the  longitudinal  cavity. 

I.  Melngailis 


B.  PROPERTIES  OF  PbSe  DIODE  LASER 
1 .  Introduction 

Semiconductor  diode  lasers  are  characterized  by  population  inversion  caused  by  minority 
carrier  injection  across  a  forward  biased  p-n  junction.  Diode  laser  action  was  first  obtained  in 
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GaAs  and  subsequently  in  other  III-V  compounds  and  alloys.  These  developments  are  reviewed 

in,  for  example.  Ref.  10.  Observation  of  diode  laser  action  in  PbTe  (Ref.  11)  broadened  the 

scope  of  laser  materials  to  include  the  lead  salts,  and  in  this  Solid  State  Research  Report  we 

12 

discuss  the  next  addition  to  the  field,  the  PbSe  diode  laser.  Laser  action  in  PbSe  occurs  at 
8.5  microns  and  extends  the  long  wavelength  limit  for  infrared  diode  lasers  beyond  the  previous 
value  of  6.5  microns  set  by  PbTe. 

It  is  significant  that  the  PbSe  laser  line  occurs  in  the  8-  to  14-micron  atmospheric  window. 

Attenuation  of  infrared  radiation  due  to  scattering  in  most  types  of  haze  resulting  from  dust, 

-4 

soot,  smoke,  etc.  (excluding  water  droplets)  varies  as  X  over  the  portion  of  the  spectrum 

13 

covered  by  diode  lasers.  Thus,  the  PbSe  laser  at  8.5  microns  may  ultimately  be  more  suitable 
for  applications  involving  transmission  through  the  atmosphere  than  the  GaAs  laser  at  less  than 

14 

one  micron.  Since  the  magnitude  of  the  energy  gap  in  PbSe  is  strongly  affected  by  temperature, 
magnetic  field/5  and  pressure ^  the  laser  can  be  tuned  over  a  significant  range. 

Since  PbSe  is  new  to  the  modern  field  of  single-crystal  devices  and  since  many  of  the  fabri¬ 
cation  techniques  for  the  PbSe  diode  lasers  differ  markedly  from  those  used  for  the  III-V  diode 
ls.sers,  the  properties  of  PbSe  and  the  details  of  laser  fabrication  will  be  extensively  discussed. 
Results  of  spectral  measurements  for  diode  currents  above  and  below  the  threshold  for  laser 
action  are  presented  and  preliminary  results  on  the  shift  in  the  spectral  position  with  magnetic 
field  are  described. 

2.  Properties  of  PbSe 

PbSe  possesses  a  rock  salt  crystal  structure,  unlike  the  more  familiar  III-V  compounds, 

which  have  the  zinc -blende  structure.  PbSe  is  able  to  exist  as  a  stable  compound  under  fairly 

wide  deviations  from  the  stoichiometric  composition  of  equal  numbers  of  lead  and  selenium  atoms. 

Furthermore,  since  every  excess  lead  atom  in  the  lattice  gives  rise  to~an  electron  donor  and 

every  excess  selenium  atom  to  an  acceptor,  the  carrier  concentration  and  type  can  be  controlled 

by  adjusting  the  Pb:Se  ratio  in  a  pure  crystal.  PbSe  may  be  doped  n-type  to  a  concentration  as 

19  -3 

high  as  approximately  5  X  10  electrons  cm  by  the  addition  of  excess  lead,  or  p-type  to  about 

19  -3  17 

10  holes  cm  by  the  addition  of  excess  selenium.  In  our  fabrication  procedures,  carrier 

type  and  concentration  were  controlled  in  this  manner. 

Since  the  conduction  band  minima  and  valence  band  maxima  for  PbSe  both  occur  on  the 

^ill>  axes  at  the  Brillouin  zone  surface  (L-points),  it  is  a  direct-gap  semiconductor  with  a 

four-fold  minimum  energy  gap  in  k-space.  This  band  structure  is  in  contrast  to  that  for  the 

III-V  compounds  which  possess  single  minimum  energy  gaps  at  the  zone  center  (T -point).  For 

PbSe,  the  surfaces  of  constant  energy  near  the  edges  of  both  bands  are  prolate  spheroids  with 

major  axes  in  the  ^111^  directions  and  anisotropy  ratios  of  approximately  2.0.  The  dependence 

o::  energy  on  k  is  nonparabolic,  hence  the  effective  masses  are  concentration  dependent.  The 

15 

transverse  mass  at  the  band  edges  for  both  electrons  and  holes  is  0.049  electron  mass  units.  ' 

Carrier  mobilities  are  approximately  equal  for  electrons  and  holes.  In  a  well-prepared 

32-1-1  42-1-1 

crystal,  the  mobility  may  be  10  cm  volt  sec  at  room  temperature,  3.5  X  10  cm  volt  sec 

at  77°K  and  5  X  105  cm2  volt-1  sec-1  at  4.2'K. 
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Fig.  1-3.  Annealing  of  single  crystal  FbSe.  Temperature  profile  along 
sealed  quartz  tube  during  annealing  is  shown.  Upper  and  lower  temper¬ 
atures  are  654°  and  195°C.  A  one-week  anneal  at  these  temperatures 
produces  n-type  PbSe  with  3X  10^  electrons  cm“^. 


VAPOR  SOURCE 


Fig.  1-4.  Formation  of  diffused  p-n  junctions  in  PbSe. 
PbSe  crystals  and  a  two-phase  PbSe,  either  lead-  or 
selenium-rich,  vapor  source  are  in  separate  sections  of 
a  sealed  quartz  ampoule. 
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3.  Annealing  and  Fabrication  Techniques 

19 

The  desired  carrier  type  and  concentration  in  PbSe  was  obtained  by  an  annealing  process 
utilizing  the  property  that  for  a  given  crystal  temperature,  the  equilibrium  Pb:Se  ratio  is  deter¬ 
mined  by  the  selenium  vapor  pressure  which,  in  turn,  is  controlled  by  the  temperature  of  a 
charge  of  pure  selenium.  As  shown  in  Fig.  1-1,  2-mm  thick  wafers  of  PbSe,  which  had  been 

previously  cut  from  a  Bridgman  grown  n-type  boule  having  an  as-grown  electron  concentration 
18  -3 

of  about  3  X  10  cm  ,  and  a  charge  of  pure  selenium  were  placed  at  opposite  ends  of  a  quartz 

tube.  The  tube  was  evacuated,  back-filled  with  argon  to  a  pressure  of  100  torrs  and  sealed  off 

before  being  placed  in  a  two-zone  furnace  with  the  temperature  profile  shown  in  the  figure.  The 

17 

parameters  given  in  the  caption  of  Fig.  1-3  resulted  in  an  n-type  crystal  with  3  X  10  electrons 

-3  2  -1  -1 

cm  and  a  77°K  Hall  mobility  of  36,000cm  volt  sec  .  For  wafer  and  selenium  temperatures 

17  -3 

of  659°  and  211°C,  respectively,  a  p-type  crystal  resulted  with  5  X  10  holes  cm  and  a  Hall 
2  -1  -1 

mobility  of  33,000cm  volt  sec  at  77°K.  Annealing  times  were  of  the  order  of  a  week.  The 

high-temperature  equilibrium  conditions  were  frozen  in  by  quenching  to  room  temperature.  Both 

n-  and  p-type  annealed  PbSe  has  been  used  for  fabricating  diode  lasers. 

20 

In  preparation  for  forming  planar  diffused  p-n  junctions,  rectangular  parallelepipeds  with 
linear  dimensions  less  than  1  mm  were  obtained  from  the  annealed  wafers  by  cleaving  on  {100} 
planes.  The  cleaved  samples  were  placed  in  one  section  of  a  quartz  ampoule  and  a  vapor  source 
in  another  section,  as  shown  in  Fig.  1-4.  To  form  an  n-layer  on  p-type  PbSe  the  vapor  source 
was  a  lead-rich,  two-phase  ingot  of  PbSe;  and  for  a  p-layer  on  n-type  material  it  was  a  two-phase, 
selenium-rich  ingot,  the  ingots  in  both  cases  being  ground  to  a  coarse  powder.  The  loaded 
ampoule  was  evacuated,  sealed  off,  heated  in  a  single-zone  furnace  and  quenched  in  air  to  room 
temperature.  Diffusion  parameters  for  several  concentrations  of  PbSe  are  given  in  Table  1-1. 

It  is  interesting  to  note  that  junctions  formed  in  this  manner  can  be  both  intrinsic  and  impurity- 
free.  The  techniques  of  impurity  doping  used  in  elemental  and  III-V  semiconductors,  on  the 
other  hand,  produce  junction  regions  that  are  intrinsic  through  compensation. 

Diodes  in  the  form  of  rectangular  parallelepipeds  were  formed  from  the  diffused  material 
by  further  cleaving  on  {l00}  planes.  Typical  dimensions  are  0.40  X  0.18  X  0.18  mm,  the  longer 
dimension  being  the  distance  between  the  parallel  reflecting  faces  which  define  the  Fabry-Perot 
structure  for  laser  action.  Ohmic  contact  to  p-type  PbSe  was  made  with  evaporated  gold,  and 
to  n-type  material  with  indium  or  an  Sn-In  alloy.  The  lasers  were  mounted  in  low-inductance, 
microwave -type  packages  as  shown  by  the  artist's  representation  in  Fig.  1-5.  Packaging  for 

cryogenic  operation  is  generally  not  difficult  since  the  thermal  expansion  coefficients  of  the  lead 
2 1 

salts  are  comparable  to  those  of  easily  worked  metals,  such  as  copper. 

22 

It  has  been  found  at  the  Laboratory  and  by  other  workers  that  lead  salt  crystals  may  be 
damaged  in  the  bulk  by  the  sawing,  grinding,  and  polishing  techniques  that  are  commonly  used 
with  other  semiconductors.  For  this  reason,  shaping  operations  in  the  fabrication  of  PbSe  lasers 
have  been  restricted  to  cleaving,  spark-cutting,  and  etching. 

4.  Spectral  Measurements 

The  diode  lasers  were  mounted  on  a  conduction-cooled  block  of  a  liquid  helium  dewar  equipped 

23 

with  a  BaF^  window.  The  block  temperature  was  monitored  with  a  GaAs  diode  thermometer. 
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TABLE  1-1 

DIFFUSION  PARAMETERS  FOR  PRODUCING  p-n  JUNCTIONS 

IN  SINGLE-CRYSTAL  PbSe  FOR  A  DIFFUSION  TIME  OF  ONE  HOUR. 

SOURCE  INGOT  IS  TWO-PHASE  PbSe  WITH  Pb:Se  ATOMIC  PERCENT  RATIO 

OF  51:49  FOR  p-TYPE  SUBSTRATE  AND  49:51  FOR  n-TYPE  SUBSTRATE. 

Crystal  No. 

Substrate  Carrier 

Type  and  Concentration 
(cm-3) 

Diffusion 

Temperature 

(°C) 

Junction  Depth 
(cm) 

29-63  C-2  (annealed) 

p  =  1 .5  X  1017 

350 

<1  x  io'3 

29-63  C-2  (annealed) 

p  =  1.5X  1017 

375 

5  X  1 0-3 

29-63  L-4  (annealed) 

p  =4.7X10,? 

530 

7  X  10-3 

1 9733-3 1  (unannea  1  ed) 

p  =  1.9X1018 

650 

9  X  1  o”3 

29-63  0-4  (annealed) 

n  =3.1  X  1017 

450 

4  X  1 0_3 

19733-32  (unannealed) 

n  =4.4X  1018 

550 

8X  10'3 
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Fig.  1-5.  (a)  Artist's  sketch  of  low-inductance  package  for  PbSe  diode  laser:  (b)  Details 
of  the  mounting  of  a  diode  laser  in  its  package.  Laser  shown  was  fabricated  from  p-type 
base  material . 
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Fig.  1-6.  Spectra  of  infrared  emission  from  a  PbSe  diode  laser 
at  12°K  with  forward  currents  below  threshold  for  stimulated 
emission.  Diode  area  =  10“^  cm^. 


WAVELENGTH  (micron*) 


CURRENT  DENSITY  (omp/cmZ) 


Fig.  1-7.  Spectrum  of  coherent  infrared  emission 
from  a  PbSe  diode  laser  at  12°K  with  forward  cur¬ 
rent  well  above  threshold  for  stimulated  emission. 
Forward  current  density  =  4,000  amp  cm-z  and 
diode  area  =  2  X  10“^  cm^. 


Fig.  1-8.  Intensity  of  infrared  emission 
from  cleaved  end  of  a  PbSe  diode  laser 
at  12°K  as  a  function  of  diode  current. 
Diode  area  =  2  X  10“^  cm^. 
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Forward  current  was  supplied  in  pulses  with  puLse  width  adjusted  to  eliminate  heating.  The  in¬ 
frared  output  from  an  end  face  of  the  diode  was  focused  with  f/4  reflecting  optics  onto  the  en¬ 
trance  slit  of  a  single-pass  monochrometer.  The  signal  was  detected  by  a  copper-activated 
germanium  detector  whose  output  was  fed  into  a  lock-in  amplifier  system. 

Figure  1-6  shows  emission  spectra  from  a  PbSe  diode  at  a  temperature  of  approximately 
12. °K  for  several  values  of  forward-current  density.  At  low  currents,  a  very  broad  curve  with 
a  peak  at  10.1  microns  is  observed.  As  the  current  is  increased,  a  narrower  second  peak  be¬ 
comes  evident  at  about  8.5  microns  and  the  amplitude  of  this  peak  increases  superlinearly  as 
the  current  is  raised.  This  superlinearly  increasing  spontaneous  radiation  at  8.5  microns  soon 
dominates  the  diode  emission,  and  above  a  threshold  current  density  laser  action  at  this  wave¬ 
length  is  observed.  The  spectrum  of  coherent  radiation  from  the  laser  at  a  current  well  above 
threshold  is  shown  in  Fig.  1-7.  The  average  wavelength  separation  between  the  modes  in  Fig.  1-7 

o 

is  287  A.  Using  the  usual  expression  for  the  axial  mode  spacing  in  the  Fabry-Perot  cavity  of 

24 

length  0.4  mm,  the  room-temperature  value  of  4.8  for  the  refractive  index  and  neglecting  the 

o 

dispersion  term,  the  mode  separation  is  calculated  to  be  200  A,  in  good  agreement  with  the 

o 

measured  value.  The  50 -A  width  of  the  main  peak  shown  in  Fig.  1-7  is  resolution  limited.  The 
absolute  intensity  of  the  coherent  radiation  in  Fig.  1-7  is,  of  course,  many  orders  of  magnitude 
larger  than  the  absolute  intensity  of  the  spontaneous  radiation  in  Fig.  1-6. 

Figure  1-8  is  a  plot  of  radiation  intensity  from  the  end  face  of  a  diode  as  a  function  of  cur¬ 
rent  density.  The  rapid  rise  at  low  currents  reflects  the  superlinear  behavior  of  the  8. 5 -micron 

_2 

peak  with  current.  The  break  in  the  curve  at  about  2,000  amp  cm  corresponds  to  the  onset  of 
laser  action.  At  higher  currents,  the  slope  becomes  linear,  indicating  that  the  quantum  ef¬ 
ficiency  reaches  a  limiting  value.  This  behavior  of  emission  with  diode  current  differs  from 
that  of  the  III-V  diode  lasers  for  which  the  spontaneous  radiation  generally  increases  linearly 
with  current,  and  superlinear  behavior  is  usually  associated  with  the  onset  of  laser  action. 

Preliminary  measurements  have  been  made  of  the  effect  of  a  magnetic  field  on  the  laser 
output.  With  the  magnetic  field  parallel  to  the  diode  current,  a  ^100^  crystallographic  direction, 

the  peak  of  the  laser  emission  in  which  cavity  mode  structure  was  not  resolved  was  observed 

—8 

tc  shift  linearly  at  the  rate  of  7.1  X  10  ev/gauss  or  17  Mcps/gauss.  Results  from  interband 

15 

rr.agneto -absorption  studies,  predict  this  magnetic  shift  for  emission  associated  with  band-to- 

band  transitions.  The  magnetic  field  also  reduced  the  threshold  current  density  from  its  zero- 

-2  -2 

field  value  of  2,000  amp  cm  to  500  amp  cm  at  a  field  of  about  lOkgauss.  For  fields  above 

this  value,  the  threshold  increased  slowly  with  increasing  field.  The  initial  decrease  in  thresh- 

25 

old  appears  similar  to  the  effect  observed  with  InSb  by  Phelan,  et  al. 

5.  Conclusions 

The  PbSe  diode  laser,  with  a  coherent  emission  peak  at  8.5  microns  at  12°K,  has  been 
developed.  Fabrication  techniques  are  based  on  controlling  carrier  type  and  concentration  by 
adjusting  the  Pb:Se  ratio,  rather  than  by  adding  impurities.  The  diode  emission  below  threshold 
exhibits  two  spectral  peaks,  one  with  a  maximum  near  8.5  microns  which  increases  superlinearly 
with  current  and  another  with  a  maximum  at  10.1  microns  which  increases  slowly  with  current. 
Laser  action  associated  with  the  8.5-micron  peak  occurs  above  a  threshold  current  density  of 
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2, 000  amp  cm  .  The  unresolved  emission  peak  shifts  linearly  to  higher  energies  with  ^100^- 

-8 

oriented  magnetic  field  intensity  at  the  rate  of  7.1  X  10  ev/gauss,  or  17  Mcps/gauss,  which  is 
the  expected  value  if  the  emission  is  associated  with  band -to -band  transitions.  The  threshold 
current  density  decreases  to  a  fraction  of  its  zero-field  value  in  a  magnetic  field  of  approximately 
lOkgauss,  then  rises  slowly  with  increasing  field.  The  PbSe  laser  should  eventually  be  of  partic¬ 
ular  interest  for  communications  since  its  emission  is  in  the  8-  to  14-micron  atmospheric  win¬ 
dow,  a  spectral  region  of  high  atmospheric  transparency  where  attenuation  due  to  scattering  by 
haze  is  generally  low.  j  p  Butler 

A.  R.  Calawa 
R.  H.  Rediker 
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II.  LASER  RESEARCH 


A.  RAMAN  LASER  PROGRAM 

1.  Stimulated  Brillouin  Emission  at  90°  to  Ruby  Beam 

Stimulated  Brillouin  laser  action  in  nitrobenzene  has  been  observed  at  90°  to  the  ruby  beam 
using  the  same  physical  arrangement  as  was  used  for  the  90°  stimulated  Raman  experiment 
described  previously.  Slight  changes  were  made  to  isolate  the  effects  from  one  another. 

o 

If  the  mirrors  used  in  the  90°  cavity  reflect  at  6940  as  well  as  7650  A,  then  both  effects  are 
observed  simultaneously,  and  in  fact,  S^  Brillouin  and  S^,  S^,  and  Raman  emission  have 
been  observed  on  the  same  1Z  infrared  plate.  However,  if  mirrors  are  used  that  reflect  6940 

o 

but  not  7650  A,  the  Raman  emission  is  not  observed.  This  enables  a  study  of  the  Brillouin 
emission  without  fogging  the  infrared  film  with  Raman  radiation.  Conversely,  if  the  mirrors 

o 

reflect  7650  and  not  6940  A,  then  the  Brillouin  emission  is  not  present. 

Study  with  an  external  etalon  crossed  with  a  prism  shows  that  the  Brillouin  emission  is 

- 1  - 1 

about  0.03  cm  or  1  kMcps  wide,  whereas  the  Raman  emission  extends  over  1.5  to  3  cm 

There  is  some  evidence  in  these  pictures  of  quenching  of  the  Brillouin  effect  by  a  strong  Raman 
effect. 

The  90°  nonstimulated  Brillouin  shift  in  nitrobenzene  was  measured  with  a  He-Ne  laser  by 

2 

the  same  method  used  by  Chiao  and  Stoicheff.  The  shift  was  5  kMcps,  as  shown  in  Fig.  II— 1 . 

The  180°  shift  calculated  from  this  is  7.1  kMcps.  This  number  agrees  closely  with  that  obtained 
in  180°  stimulated  Brillouin  emission  with  ruby  by  Garmire.3 

External  etalon  photographs  of  the  90°  stimulated  Brillouin  output  (Fig.  II-2)  show  one 
Brillouin  mode  designated  B^  and  sometimes  another  ring,  designated  B^,  between  6.8  and 
7.1  kMcps  away  from  B^.  None  of  the  primary  ruby  radiation  appears  in  Figs.  II- 1  and  II-2. 

In  order  to  see  the  relationship  between  the  ruby  modes  and  the  first  Brillouin  (B^)  ring, 
some  of  the  0°  ruby  energy  was  fed  around  by  a  mirror  and  added  to  the  90°  emission  analyzed 
by  the  external  etalon.  A  typical  picture  is  shown  in  Fig.  II-3.  A  study  of  many  of  these  pic¬ 
tures  shows  that,  in  general,  no  one  ruby  mode  is  exactly  5  kMcps  away  from  the  first  Brillouin 
ring.  Therefore,  one  of  the  ruby  modes  alone  could  not  have  produced  the  B^  ring  by  itself. 

Each  ruby  mode  produces  a  Lorentzian-shaped  gain  line  spaced  5  kMcps  from  itself.  The 
composite  gain  line  is  a  sum  of  these,  and  its  shape  depends  on  the  relation  between  the  spacing 
of  the  ruby  modes  and  the  bandwidth  of  the  BriLouin  scattering  process.  An  upper  limit  to  the 
w..dth  of  the  process  is  known  from  the  fact  that  the  nonstimulated  emission  is  no  wider  than  the 
0.6-kMcps  gas  laser  that  produced  it.  The  width  of  a  single  ruby  mode  is  also  approximately 
0.6  kMcps  so,  in  this  case,  the  Lorentzian  gain  line  will  be  of  this  width,  and  the  composite 
line  will  have  peaks  which  are  double  the  amplitude  of  the  valleys. 

A  study  of  many  far-field  patterns  indicates  that  the  system  usually,  but  not  always,  oscil¬ 
lates  in  an  on-axis  mode  of  the  cavity  with  radiation  extending  about  1  kMcps  around  the  center 
frequency.  Hence,  the  resonant  cavity  of  4.7-kMcps  interorder  spacing  usually  picks  a  fre¬ 
quency  on  the  composite  gain  line  at  which  to  oscillate,  which  is  not  always  a  peak  of  the  line. 
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I  -82  -  571 1] 


Fig.  1 1  —  1 .  90°  incoherent  Brillouin  scatter 
using  a  1 3-mw  6328  A  He-Ne  laser  source. 
Etalon  spacing  is  0.25  inch.  Center  ring  of 
each  group  is  laser  light.  Outer  ring  is  anti- 
Stokes,  inner  ring  is  Stokes  scatter.  Magni¬ 
tude  of  the  shift  is  5  kMcps. 


Fig.  11-2.  External  etalon  analysis  of  90° 
Brillouin  output.  Etalon  spacing  is  0.  25 
inch. 


|  3-82-3715-1] 
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R's  ARE  RUBY  MODES  |  3-82*3716-1  1 

B,  IS  90*  BRILLOUIN  EMISSION  1 - L 


G's  DENOTE  POSITIONS  OF  RUBY 
MODES  SHIFTED  BY  5kMcps 


r 


GAIN  LINE 
PRODUCED  BY 
90*  BRILLOUIN 
SHIFT 


Fig.  11-3.  Relationship  of  ruby  frequency  modes 
to  90°  Brillouin  output. 
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Occasionally,  the  far -field  pattern  is  dark  in  the  center  indicating  off-axis  oscillation  in  the 

main  B  .  line. 

*  4 

Bloembergen  and  Shen  have  shown  that  different  ruby  pump  modes  couple  to  each  other 
through  the  Raman  effect  and  produce  cross  terms  in  the  Raman  gain  expression  that  result  in 
more  gain  than  if  the  same  energy  were  put  into  a  single  ruby  mode.  This  same  argument, 
applied  to  the  Brillouin  case,  might  explain  the  difficulty  in  doing  the  90°  Brillouin  experiment 
with  a  single -frequency  ruby  pump.  Threshold  was  never  reached  with  a  single  mode  of  0.2  Mw 
using  a  passive  Q-switched  ruby.  It  is  estimated  that  the  rotating  mirror  system  gives  about 
0.5  Mw  in  its  stronger  modes  so  that  it  is  impossible  to  rule  out  that  more  power  in  a  single 
mode  might  reach  threshold.  All  attempts  to  produce  a  single  mode  from  the  rotating  mirror 
system  failed. 

This  experimental  configuration  enables  the  gain  of  the  Brillouin  scattering  medium  to  be 

measured  from  the  conditions  at  threshold.  Reamplification  of  the  scattered  light  in  the  ruby, 

which  occurs  in  the  180°  backscattered  experiments,  does  not  occur  in  this  90°  configuration. 

The  threshold  for  this  stimulated  Brillouin  effect  is  about  3  Mw  from  our  ^-inch  rod.  The  rod 

is  judged  to  be  lasing  strongly  over  5  mm  of  its  diameter.  This  is  the  "path  length"  for  gain 

across  the  cylindrical  focus.  The  cylindrical  lens  concentrates  the  beam  10:1  in  one  dimension 

2 

giving  an  active  area  of  0.05  cm  .  The  gain  at  threshold  is  estimated  to  be 

a  =  *  .  ^  n  5* -99*  -  o.02cm.  *  at  60  Mw/cm2  incident  flux 
l  O.b  cm 

At  5-Mw  ruby  power,  or  100-Mw/cm  power  density,  the  90°  Brillouin  output  is  estimated 
to  be  1  kw.  Inside  the  cavity  of  99.1-percent  mirrors,  the  power  is  then  about  100  kw.  The 
threshold  for  this  90°  experiment  is  somewhat  lower  than  that  reported  for  the  0°  experiment 
us:mg  a  spherical  focus,  but  no  resonator'*  around  the  active  medium. 

There  is  reason  to  wonder  how  a  100 -kw  (inside  the  cavity)  B^  line  can  produce  a  B^  line 
by  the  process  of  180°  Brillouin  scattering  when  the  threshold  for  the  production  of  the  B^  line 
is  3  Mw  or  an  average  of  about  |Mw  per  (strong)  mode  of  the  ruby.  One  possible  explanation 
is  that  the  threshold  for  a  180°  backward  wave  process  is  lower  than  for  the  90°  process,  or 
that  there  is  some  additional  gain  from  additional  ruby  modes  spaced  5  kMcps  above  the  B^ 
line.  The  B^  line  is  shown  in  Figs.II-2  and  II-4. 


Fig.  11-4.  External  etalon  analysis  of  sum  of  90° 
Brillouin  output  and  a  small  amount  of  ruby  pump 
light  in  a  case  when  a  second  Brillouin  line  is 
present. 


R's  ARE  RUBY  MODES 
B's  ARE  90*  BRILLOUIN  EMISSION 
6  s  DENOTE  RUBY  MODES 
SHIFTED  BY  5kMcps 
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The  ability  of  the  cavity  of  4.7-kMcps  interorder  spacing  to  support  two  modes,  and  B^, 
7  kMcps  apart,  can  be  explained  by  the  fact  that  a  plane  parallel  cavity  with  an  aperture  as 
large  as  its  spacing  can  have  walkoff  modes  of  considerable  angle.  If  the  B^  mode  is  assumed 
to  oscillate  on  the  axis  of  the  cavity,  there  will  be  well  over  100  passes  before  the  B^  mode, 
which  is  7.0  —  4.7  =  2.3  kMcps  from  an  on-axis  resonance,  walks  off  the  mirrors  of  2-cm 
aperture. 

5  1 

The  stimulated  Brillouin  effect  might  be  even  better  for  multiple  pumping  schemes  '  to 
obtain  high  power,  with  high  spatial  mode  purity  than  is  the  Raman  effect.  The  phonon  created 
is  of  much  lower  energy  and  hence  results  in  much  less  local  heating  of  the  medium  to  cause 
index  of  refraction  changes  and  a  degradation  of  optical  quality. 


2.  Stimulated  Raman  Emission  at  90°  to  Ruby  Beam 

Some  additional  work  has  been  done  on  the  90°  Raman  experiment.  The  90°  Raman  output 

in  nitrobenzene  was  analyzed  with  an  external  etalon,  and  no  structure  was  observed  in  the 

_  \ 

rings.  The  output  was  a  continuum  of  about  2-cm  width.  It  was  conjectured  that  the  plane 
parallel  Raman  cavity  was  oscillating  in  walkoff  modes  which  could  then  range  in  frequency 
continuously  around  the  center  frequency. 

To  show  this,  the  far  field  pattern  of  the  90°  Raman  output  was  focused  in  a  plane.  Then, 

a  slit  was  used  to  block  off  most  of  the  off-axis  energy.  The  remaining  energy,  when  fed  into 

an  etalon,  showed  a  large  number  of  modes  separated  by  the  interorder  spacing  of  the  Raman 

_  \ 

cavity  covering  a  range  of  approximately  2  cm  .  Jane  H  Dennis 


3.  Stimulated  Stokes  and  Anti-Stokes  Radiation  from  a  Point  Source 

We  have  considered  a  model  for  the  Stokes  and  anti-Stokes  radiation  arising  from  a  point 
source  at  the  Stokes  frequency  in  a  Raman-active  medium.  The  coupled  Fourier  transform  equa¬ 
tions  for  the  Stokes  and  anti-Stokes  fields  are  used  to  obtain  an  exact  Fourier  integral  expres¬ 
sion  for  the  Stokes  and  anti-Stokes  fields.  These  Fourier  integral  expressions  are  evaluated 
asymptotically  using  a  saddle -point  method.  It  is  found  that,  to  lowest  order  in  the  Raman 
gain  and  coupling  constants,  the  Stokes  wave  grows  in  a  spherical  wave  from  the  point  source 
with  the  same  growth  constant  in  the  phase -matched  direction  as  in  other  directions.  The  re¬ 
action  of  the  anti-Stokes  radiation  on  the  Stokes  radiation  in  the  phase-matched  direction  pro¬ 
duces  only  a  relative  suppression  in  higher  order;  the  growth,  however,  remains  exponential. 

The  anti-Stokes  radiating  wave  in  the  Raman-active  medium  is  exponentially  damped,  be¬ 
cause  the  source  of  phase -matched  anti-Stokes  radiation  is  a  narrow  cone  of  polarization  at  the 
anti-Stokes  frequency  associated  with  the  phase -matched  Stokes  radiation,  and  as  the  anti-Stokes 
radiation  propagates  away  from  this  active  cone,  it  is  damped  by  the  anti-Stokes  damping  of  the 
active  medium.  Thus,  the  only  anti-Stokes  radiation  observed  externally  is  the  radiation  pro¬ 
duced  by  the  polarization  at  the  anti-Stokes  frequency  within  the  phase -matched  Stokes  cone  and 
a  boundary  of  the  Raman-active  medium.  ^  ^  Haus* 

p!  l!  Kelley 
H.  J.  Zeiger 


*  Department  of  Electrical  Engineering,  M.I.T. 
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B.  CW  TWO-PHOTON  PHOTOELECTRIC  EFFECT 


We  have  observed  continuous  double -photon  photoelectric  emission  from  Cs~Sb  using  the 

^  6 

1.153-micron  radiation  from  a  He-Ne  gas  laser.  The  first  observation  of  this  effect  utilized 
a  pulsed  solid-state  laser.  In  our  experiment,  the  laser  light  was  focused  upon  the  cathode  of 
an  RCA  1P21  photomultiplier  tube  by  a  4.8-cm  focal-length  lens.  With  the  laser  output  set  at 
1.35  mw,  calibrated  attenuators  were  employed  to  vary  the  power  incident  upon  the  cathode  sur¬ 
face.  Reflection  losses  between  the  exit  mirror  of  the  laser  and  the  cathode  amounted  to 
30  percent. 

Figure  II-5  shows  preliminary  results  obtained  with  the  laser  oscillating  in  a  higher -order 

mode  whose  pattern  occupied  an  area  approximately  5  mm  in  diameter  at  the  plane  of  the  lens. 

In  computing  the  current  and  power  densities  of  Fig.  II- 5  it  was  assumed  that  the  focal  spot  had 
-5  2 

an  area  of  10  cm  ,  which  is  approximately  five  times  the  theoretical  diffraction-limited  focal 
spot  area  of  the  TEM^  mode.  Unfortunately,  the  focal  spot  was  small  enough  to  reveal  non¬ 
uniformities  across  the  cathode  surface,  thus  hampering  reproducibility  of  the  data.  Also,  a 
fraction  of  the  convergent  laser  beam  was  intercepted  and  scattered  by  the  grill  wires  used  for 
electrostatic  focusing  of  the  emitted  electrons. 

The  quadratic  portion  of  the  curve  in  Fig.  II-5  represents  two-photon-induced  emission  as 
discussed  by  Sonnenberg,  et  al.^  The  linear  portion  is  characteristic  of  a  one -photon  process, 

7 

ancl  similar  behavior  has  been  seen  by  Teich,  et  al.,  in  a  two-photon  experiment  with  sodium 
metal.  In  the  present  experiment,  the  linear  portion  may  be  attributed  to  the  finite  occupancy 
of  the  conduction  band  of  Cs^Sb  and  to  the  presence  of  occupied  impurity  levels  and  other  states 
in  the  forbidden  band  lying  within  about  0.7  ev  of  the  bottom  of  the  conduction  band. 

R.  A.  Soref 


Fig.  11-5.  Photoemissive  response  of  Cs^Sb 
cathode  to  CW  1.  153-micron  radiation. 
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|-SS  -3836  1 


Fig.  11-6.  Photoconductive  response  of  semiconductor  samples  to  unfocused  6943-A  laser  radiation: 
(a)  High-purity  CdS,  pQ  =  5  X  10?  ohm-cm;  (b)  CdS,  pQ  =  3  X  10^  ohm-cm;  (c)  CdS,  =  5  ohm-cm; 
(d)  High-purity  ZnS,  po  =  3  X  10^  ohm-cm;  (e)  High-purity  ZnSe  polycrystal ,  po  =  2  X  10^  ohm- 
cm.  Upper  trace:  laser  intensity.  Lower  trace:  photocurrent,  time  scale  =  50- psec  div“^ . 
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C.  TWO-PHOTON  PHOTOCONDUCTIVITY 

If  a  semiconductor  is  illuminated  with  intense  optical  radiation  of  frequency  double  - 

-1 

phDton  absorption  will  occur  provided  a;  <  h  <  2a>,  and  the  hole-electron  pairs  thereby  cre¬ 
ated  will  alter  the  electrical  conductivity  of  the  crystal.  We  have  attempted  to  observe  " second- 
order"  photoconduction  of  this  sort  in  monocrystals  of  CdS,  ZnS,  and  ZnSe  using  ruby  laser 

radiation.  Representative  preliminary  results  are  shown  in  Fig.II-6.  The  largest  laser  spikes 

22  -2  -1 

in  Fig.II-6  correspond  to  a  peak  flux  of  about  10  photons  cm  sec  incident  upon  the  samples. 
Tr  apping  effects  may  have  caused  the  steady  DC  rise  in  photocurrent  over  a  200-psec  interval 
in  samples  (a)  and  (b),  although  this  behavior  is  not  fully  understood  at  present.  Dielectric  re¬ 
laxation  may  be  ruled  out  on  the  basis  of  the  low  resistivity  of  sample  (b). 

As  the  output  energy  of  the  laser  was  increased  from  0.1  to  0.5  joule  per  flash,  it  was  found 
that  the  time-integrated  photocurrent  of  sample  (a)  increased  as  the  cube  of  the  time-integrated 
laser  intensity,  whereas,  for  samples  (d)  and  (e),  the  integrated  photocurrent  increased  linearly 
with  the  integrated  laser  intensity.  The  cubic  response  has  not  yet  been  explained.  The  linear 
response  could  be  attributed  to  two-photon  absorption,  provided  that  bimolecular  recombination 
across  the  gap  is  the  dominant  recombination  mechanism.  However,  single-photon  impurity 
absorption  can  also  yield  a  linear  response  if  a  linear  recombination  process  dominates.  There¬ 
fore,  additional  data  are  required  to  decide  how  many  charge  carriers  (if  any)  were  produced 
via  two-photon  absorption.  ^  ^  Soref 

M.  M.  Litvak 


D.  GAS  LASER  PROGRAM 

1.  Measurement  of  Amplitude  Noise  in  Optical  Cavity  Masers 

We  report  here  measurements  on  the  amplitude  noise  of  a  cavity-type  gaseous  optical 

o 

maser  oscillator  operating  at  6328  A  above  and  below  threshold.  We  compare  our  results  above 

g 

threshold  with  the  theory  of  noise  in  van  der  Pol  oscillators,  which  is  the  proper  one  to  use  for 
q 

the  maser  oscillator,  and  the  measurement  below  threshold  with  the  linear  theory  of  maser 
10 

amplifiers.  Because  of  the  quantitative  agreement  we  conclude  that  the  noise  is  due  to 
spontaneous  emission. 

We  use  the  equivalent  circuit  of  Fig.  II- 7 .  The  conductance  G  =  cjoC/Q  represents  the 
resonator  loading  by  the  outside  space. 

is  the  (usually  negative)  conductance  of  the  maser  material.  Above  threshold,  the  magnitude 
of  the  negative  conductance  Gm  decreases  quadratically  with  voltage.  The  uniform  spectral 
density 

iZ/Av  =  4Gmhy  (n2/g2)/[(n1/gj)  -(n2/g2)] 

of  the  noise  current  generator  is  chosen  so  that  the  power  emitted  by  the  cavity  below  threshold 

10 

satisfies  the  well-known  linear  theory  (n^  and  n^  are  the  populations  of  the  lower  and  upper 
levels,  and  g^  and  g^  are  their  degeneracies).  Above  threshold,  the  dependence  of  the  noise 
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Fig.  11-7.  Equivalent  circuit  of  van  der  Pol 
oscillator  with  noise. 


power  upon  the  population  inversion  is  assumed  to  be  the  same  as  that  below  threshold.  Using 
the  noise  power  predicted  from  the  equivalent  circuit,  and  the  photoemission  theory  of  Kelley 
and  Kleiner,  one  obtains  for  the  spectrum  of  the  anode  current  in  a  photomultiplier  illuminated 
by  the  maser,  operating  in  a  single  mode  near  threshold,  both  above  and  below. 


S(i>)  =  2Aen  +  4AeI  na 
a  a 


1  +  <£>2 


S  +  S 


(1) 


In  Eq.  (1),  A  is  the  photomultiplier  gain,  e  is  the  electron  charge,  1^  is  the  photomultiplier 
anode  current,  T  is  the  shot -noise  enhancement  factor  caused  by  secondary  emission,  77  is  a 
loss  factor  including  the  quantum  efficiency  and  all  other  losses  between  the  maser  output  and 
photomultiplier,  or  =  (n2/g2)/[  (n2/g2)  -  (n^g^J  ,  and 


Av 


(2) 


with  i/q,  the  frequency  of  the  mode.  Note,  however,  that  below  threshold  contributions  of  other 
modes  than  the  most  strongly  excited  axial  mode  become  more  important  as  the  excitation  is 
decreased;  Sg  and  denote  the  shot -noise  and  excess  noise,  respectively. 

The  rate  n  of  photon  emission  from  the  maser  above  threshold  gives  information  on  the 
coefficient  of  the  nonlinear  part  of  the  conductance 


n 


(G^  +  G)  g 
2aG°  hy 


Below  threshold,  one  finds 


n 


a 


2ir(Avo)Z 


Av 


(3) 


(4) 


In  the  experiments,  an  internal  mirror  laser  of  near  hemispherical  geometry,  with  a  mirror 
spacing  of  50  cm,  was  used.  The  inside  diameter  of  the  quartz  tube  was  3.5  mm.  The  maser  was 
DC-excited  (gas  discharge  fluctuations  were  suppressed  by  optimized  external  circuitry  and  a 
magnet  near  the  discharge),  and  was  stabilized  to  less  than  1  percent  long-term  drift  in  intensity 
by  a  feedback  circuit  with  a  time  constant  of  i/50  sec,  by  using  the  light  from  one  of  the  maser 
mirrors.  The  light  through  the  other  mirror  was  fed  partly  to  one  photomultiplier  connected  to 
a  counter,  and  partly  to  another  connected  to  a  spectrum  analyzer.  The  spectrum  was  studied  in 
the  0-  to  17-Mcps  range  Different  spectrum  analyzers  were  used  for  the  ranges  of  the  frequency 
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spectrum.  Variation  in  the  scanning  velocities  was  from  0  to  4  cps/sec  at  the  low-frequency 
end  to  0  to  1000  cps/sec  in  the  high-frequency  range  of  the  spectrum.  The  corresponding  IF 
bandwidths  ranged  from  a  minimum  of  10  cps  to  a  maximum  of  10,  000  cps.  The  spectrum 
analyzer  output  was  fed  into  a  spectral  density  analyzer  with  integration  times  of  1  or  5  sec, 
and  then  to  a  chart  recorder. 

The  composite  experimental  results  are  shown  in  Fig.  II-8(a-b).  The  curves  are  labeled 

—  -i 

by  their  corresponding  photoelectron  rate  N  (in  sec  )  at  the  photocathode.  Above  threshold, 

7  of  the  12  dynodes  of  the  photomultiplier  were  shorted  to  the  anode.  Lorentzian  curves  can  be 
fitted  with  very  good  accuracy  onto  the  curves  above  threshold  and  to  the  top  curves  below  thresh¬ 
old.  The  two  curves  with  the  highest  ratios  Se/Sg  and  with  the  narrowest  bandwidths  represented 
the  highest  power  below  threshold  and  the  lowest  power  above  threshold  at  which  our  maser 
could  operate,  stably  and  reproducibly,  for  several  months.  The  small  spike  modulation  at 
approximately  l68kcps,  both  above  and  below  threshold,  is  due  to  residual  gas  plasma  oscil¬ 
lation. 

The  experimentally  observed  spectral  curves  of  the  photomultiplier  current,  predicted 
quantitatively  by  Eq.  (1)  and  (2),  show: 

(a)  The  ratio  Se/Ss  (near  zero  frequency),  is  proportional  to  the  maser  power 
P  below  threshold,  and  inversely  proportional  to  P  above  threshold.  De¬ 
viation  from  the  l/P  dependence  occurs  only  farther  below  threshold  when 
more  than  one  linearly  polarized  mode  contributes  to  the  output. 

(b)  The  bandwidth  of  the  excess  noise  is  inversely  proportional  to  P  below 
threshold  and  directly  proportional  to  P  above  threshold. 

(c)  The  measured  ratio  Se/Ss  from  the  N  =  3.68  x  10^  curve  above  threshold, 
yields  (Ref.  13)  a  =  1.97,  using  Eq.  (1),  the  measured  values  of  17  =  4.9  X 
10"3,  r  =  1.3,  and  taking  Ai/q  =  4.7  X  105cps  as  computed  from  the  trans¬ 
mittance  of  the  mirrors. 

(d)  From  Eq.  (1),  the  ratio  Se/Ss  above  and  below  threshold  is  the  same  for 
curves  of  equal  bandwidth  if,  below  threshold,  one  linearly  polarized  mode 
predominates  (as  is  the  case  in  the  operating  condition  of  the  N  =  1.64  x 

10 7  curve).  Note  that  G&  is  changed  only  by  a  fraction  o£a  percent  as 
the  laser  is  brought  from  below  threshold  operation  (the  N-=  1.64  x  107 
curve)  to  oscillating  condition  (the  N  =  3.68  X  10^  curve),  and  therefore 
a  should  not  change  perceptibly.  The  two  top  curves  in  Fig.  II-8(a-b) 
have  comparable  bandwidths.  If  one  predicts  from  the  above -threshold 
curves  the  ratio  Se/Ss  below  threshold,  taking  into  account  the  difference 
in  bandwidths,  one  finds  Se/Ss  =  4600  which  should  be  compared  with  the 
measured  value  of  4400  for  the  case  of  N  =  1.64  x  107. 

(e)  From  Eq.  (4),  one  may  predict  the  photon  emission  rate  from  the  maser 
(with  ct  =  1.97  for  the  curve  N  =  1.64  X  10 7)  to  be  n  =  3.22  x  10^  sec"1. 

The  measured  rate  was  N/rj  =  n  =  3.34  X  10^  sec"1. 

Our  measurements  of  S  /S  ,  bandwidth,  and  power  agree  with  the  theoretically  predicted 
values  within  5  to  10  percent  which  is  about  the  accuracy  of  our  absolute  power  measurements 
and  the  gain  stability  of  photomultipliers  in  a  normal  laboratory  environment. 

We  acknowledge  many  helpful  discussions  with  M.M.  Litvak,  A.  Javan,  and  A.  Szoke. 

C.  Freed 
H.  A.  Haus 
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(a)  Below  threshold  operation. 


1  10  10*  io5  I04  \£  I  <f 

FREQUENCY  (cp») 

(b)  Above  threshold  operation. 


Fig.  11-8.  Observed  photomultiplier  current  spectra. 
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2.  High-Pressure  Mercury  Laser 

The  possibility  of  a  high-pressure  mercury  arc  as  a  laser  medium  has  been  discussed  in 

3*C 

the  last  Solid  State  Research  Report.  Since  then,  it  has  been  shown  that  the  impurity  content 
in  a  discharge  can  be  controlled  well  enough  so  that  the  mercury  molecule  upper  laser  state 
can  be  pumped  to  an  appreciable  density.  Stabilization  of  the  pulsed  arc,  i.e.,  keeping  the 
constricted  arc  straight  and  parallel  to  the  tube  axis,  was  accomplished  by  maintaining  a 
mitliampere  DC  discharge  through  the  tube.  This  worked  in  a  mixture  of  argon  and  mercury 
up  to  approximately  lOOtorrs  of  mercury  and  400 -amp  pulsed  current,  at  which  point  the  arc 
again  became  unstable.  The  upper  level  population,  according  to  the  spontaneous  emission, 
had  a  greater  than  linear  growth  with  pressure  and  arc  current  up  to  the  unstable  region  as 
shown  in  Fig.  II— 9 . 


Fig.  11-9.  Power  output  In  mercury  vapor  green  band 
as  a  function  of  pulsed  arc  current  for  specified  mer¬ 
cury  partial  pressures.  Argon  partial  pressure  was 
25  torrs. 
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Estimates  of  the  gain  from  spontaneous  emission  from  the  laser  state  for  presently  achiev¬ 
able  pressure  and  current  conditions  indicate  that  more  than  0.3 -percent  loss  in  the  optical 
cavity  will  not  allow  oscillation.  Therefore,  an  optical  cavity  is  being  constructed  with  mirrors 
having  less  than  0.2  percent  estimated  loss  and  with  a  laser  tube  having  Brewster  angle  windows 
of  less  than  0.05 -percent  loss.  Special  cleaning  and  assembly  techniques  are  being  used  to 
maintain  high  quality  for  operation  at  about  250 °C. 

R.  J.  Carbone 
M.  M.  Litvak 

*  Solid  State  Research  Report,  Lincoln  Laboratory,  M.I.T.  (1964:3). 
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III.  MATERIALS  RESEARCH 


A.  CZOCHRALSKI  GROWTH  OF  SAPPHIRE  CRYSTALS 
IN  A  RESISTANCE  FURNACE 

Because  of  its  high  melting  point  (2050°C)  and  corrosiveness  toward  many  crucible  materials, 

sapphire  (Al^O^)  has  classically  been  grown  by  the  Verneuil  flame  fusion  technique.  Growth  of 

ruby  (Cr-doped  Al^O^)  crystals  by  the  Czochralski  method  has  recently  been  reported  by  Pala- 

diro  and  Roiter,*  who  used  an  induction-heated  tungsten  crucible  in  a  nitrogen  atmosphere,  and 

by  the  Linde  Company,2  which  has  not  described  the  details  of  its  procedure.  We  have  grown 

3 

sapphire  crystals  in  a  resistance  furnace  recently  developed  for  operation  in  oxidizing,  neutral, 
or  reducing  atmospheres  at  temperatures  up  to  2400 °C.  Sapphire  was  selected  as  a  trial  mate¬ 
rial  in  order  to  test  the  high-temperature  capability  of  the  furnace.  The  furnace  was  modified 
for  Czochralski  growth  by  adding  mechanisms  for  raising  and  lowering  the  seed  and  for  rotating 
the  seed  at  0  to  200  rpm.  Sight  ports  were  placed  in  the  furnace  lid  to  allow  projection  of  the 
image  of  the  melt  surface  on  a  screen  for  monitoring  growth.  The  melt  was  contained  in  an 
iridium  or  molybdenum  crucible  under  an  argon  atmosphere.  Since  the  molten  sapphire  did  not 
evaporate  appreciably,  it  was  possible  to  remove  the  zirconia  tube  used  to  protect  the  tantalum 
heater  from  oxygen  and  other  reactive  gases. 

A  spoke-like  pattern  of  light  and  dark  on  the  surface  of  the  melt  is  observed  on  the  monitor 
screen  before  the  seed  is  dipped.  This  pattern  is  probably  caused  by  convection  cells  resulting 
from  high  temperature  gradients  in  the  melt.  The  existence  of  such  gradients  is  established  by 
optical  pyrometer  measurements  which  show  that  the  temperature  of  the  crucible  wall  during 
growth  is  about  150°C  above  the  melting  point  of  sapphire.  After  the  rotating  seed  is  dipped 
into  the  melt,  the  radial  arms  of  the  spoke-like  pattern  become  skewed  due  to  the  influence  of 
the  seed  rotation  on  the  melt.  The  pattern  provides  a  good  view  of  the  growing  meniscus.  In 
general,  crystal  growth  starts  easily  but  becomes  more  difficult  to  maintain  as  the  crystal  be¬ 
comes  longer,  probably  because  the  temperature  gradient  in  the  solid  is  not  sufficient  to  stabi¬ 
lize  the  growing  interface. 

Crystals  grown  from  iridium  crucibles  are  water-white.  Two  of  these  crystals  are  shown 
in  Fig.  Ill—  1 .  Crystals  grown  from  molybdenum  crucibles  have  a  light-gray  tinge.  This  colora¬ 
tion  is  probably  not  due  to  the  presence  of  molybdenum  as  an  impurity,  since  no  molybdenum 
was  detected  in  a  mass  spectrographic  analysis  with  a  detection  limit  of  40  ppm.  Possibly  the 
coloration  is  due  to  an  oxygen  deficiency  resulting  from  reduction  by  the  molybdenum,  which 
has  a  greater  reducing  power  than  iridium.  ^  g  Reed 

R.  E.  Fahey 


B.  FORCED  CONVECTION  GROWTH  OF  IODINE  CRYSTALS 

The  results  of  initial  experiments  on  the  vapor  growth  of  iodine  crystals  by  forced  convec- 

4  5 

tion  were  reported  previously.  '  These  experiments  have  been  extended  in  order  to  obtain  addi¬ 
tional  data  concerning  constitutional  supercooling  of  the  vapor  and  the  dependence  of  the  growth 


23 


J 


Section  III 


Fig.  1 1 1-1.  Sapphire  crystals  grown  by  Czochralski  method  in  a  resistance  furnace. 


O  50  100  150 

6  re /cm) 

Fig.  111-2.  Growth  velocity  v  as  a  function  of  thermal  gradient  G  for  iodine  crystals 
grown  from  vapor  by  forced  convection.  For  a  fixed  value  of  G,  different  values  of  v 
were  obtained  by  varying  helium  flow  rate  and  difference  between  source  temperature  T0 
and  interface  temperature  Tj.  Whether  growth  interface  was  smooth  or  rough  is  indicated 
for  each  experimental  point.  Theoretical  curve  gives  values  calculated  from  Eq.(l)for 
maximum  velocity  vmQx  at  which  growth  can  occur  without  constitutional  supercooling. 
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rate  on  experimental  conditions.  The  experimental  procedure  was  essentially  the  same  as 
described  previously:  helium  gas  was  saturated  with  iodine  vapor  at  temperature  Tq,  heated 
to  a  higher  temperature  to  produce  unsaturation,  and  then  passed  around  a  glass  cold  finger 
where  condensation  occurred  at  an  interface  temperature  T^.  The  growth  rate  v  was  measured 
as  a  function  of  (Tq  —  T^)  and  the  thermal  gradient  G  in  the  gas  phase.  In  the  present  experi¬ 
ments,  T  was  90.0 °C  and  measurements  were  made  at  two  different  helium  flow  rates,  which 
'  o 

resulted  in  iodine  vaporization  rates  of  3.3  and  10.7  gram/hour. 

In  the  earlier  experiments,  it  was  found  that  under  given  growth  conditions  the  crystal- 
vapor  interface  remained  smooth  only  for  growth  velocities  below  a  certain  upper  limit.  The 
roughening  of  the  interface  which  was  observed  at  higher  values  of  v  was  attributed  to  constitu¬ 
tional  supercooling  of  the  vapor,  which  theoretical  analysis  showed  to  occur  when  v  exceeds  a 
maximum  value  (vmax)  given  by 
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pRT2 
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wnere  M  is  the  molecular  weight  of  the  condensing  species,  D  is  the  diffusion  constant  of  the 
species  in  the  gas  phase,  pQ  =  A  exp  [—  B/T]  is  the  saturated  vapor  pressure  at  absolute  tempera¬ 
ture  T,  and  p  is  the  crystal  density.  The  data  obtained  in  the  present  experiments  provide  ad¬ 
ditional  confirmation  of  this  explanation.  These  results  are  shown  in  Fig.  Ill  —  2,  where  values 
of  v  measured  under  various  growth  conditions  are  plotted  against  G.  The  theoretical  curve  in 
Fig.  Ill— 2,  which  gives  values  of  vmax  calculated  from  Eq.  (f),  is  seen  to  be  in  good  agreement 
with  the  experimental  boundary  between  smooth  and  rough  growth. 

Values  of  v  for  smooth  growth  as  a  function  of  (T  —  T^)  are  shown  in  Fig.  Ill— 3 .  Within  the 
limits  of  experimental  error,  v  is  independent  of  G  for  G  between  13°  and  140°c/cm.  The 
value  of  v  is  equal  to  the  rate  at  which  iodine  diffuses  to  the  growth  interface.  If  it  is  assumed 
that  the  diffusion  layer  and  thermal  boundary  layer  have  the  same  thickness,  it  follows  from 
elementary  boundary  layer  theory ^  that 
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where  K  is  a  dimensionless  constant  determined  by  the  geometry  of  the  system,  V  is  the  free 
si  ream  velocity  of  the  gas,  v  is  the  kinematic  viscosity  of  the  gas,  f  is  a  characteristic  length, 
and  p  °  and  p^  are  the  saturated  vapor  pressures  at  Tq  and  T.,  respectively.  The  theoretical 
curves  shown  in  Fig.  Ill— 3,  which  give  satisfactory  agreement  with  the  experimental  data,  were 
calculated  according  to  Eq.  (2)  by  adopting  4.9  as  the  value  of  K. 

According  to  Eq.  (2),  v  =  0  for  p^  =  (T^/T^)  p^  rather  than  for  p^  =  p°,  because  the  growth 
rate  is  determined  by  the  concentration  gradient  in  the  gas  phase,  which  depends  on  both  the 
pressure  and  temperature  gradients.  For  a  sufficiently  large  temperature  gradient,  Eq.  (2)  pre¬ 
dicts  that  growth  should  stop  when  T^  is  still  several  degrees  below  Tq.  Although  some  evidence 
for  this  result  was  obtained  in  the  present  experiments,  it  could  not  be  definitely  established  be¬ 
cause  of  the  uncertainty  in  measuring  T^  in  the  presence  of  large  temperature  gradients. 

The  highest  growth  velocity  achieved  in  the  present  experiments  was  4  mm/hour.  The 
velocity  was  limited  by  the  cooling  capacity  of  the  cold  finger  and  the  maximum  rate  at  which 
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Fig.  1 1 1—3-  Data  of  Fig.  Ill  —2  for  smooth  growth  replotted  to  show  growth  velocity  v 
as  a  function  of  (TQ  —  Tj)  for  iodine  crystals  grown  from  the  vapor  by  forced  convection. 
Experimental  points  do  not  indicate  values  of  G  at  which  they  were  obtained.  Theo¬ 
retical  curve  was  calculated  from  Eq.  (2). 


Fig.  111-4.  Seebeck  coefficient  a  vs  absolute  temperature  T  for  a  single  crystal  of 
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icdine  could  be  vaporized  without  channeling  the  crystal  bed  used  as  the  source  of  vapor.  Higher 
growth  rates  would  be  possible  if  the  apparatus  were  modified  to  remove  these  limitations. 

T.  B.  Reed 
W.  J.  LaFleur 


C.  CRYSTAL  GROWTH  AND  ELECTRICAL  PROPERTIES 
OF  TITANIUM  SUBOXIDES 

3  3 

Single  crystals  of  Ti^O^  (~1  cm  )  and  TiO  (1  mm  )  have  been  grown  by  the  Czochralski 
method  in  the  resistance  furnace  described  in  Sec.III-A.  The  compounds  were  prepared  by  re¬ 
ducing  well-dried  TiO^  powder  with  appropriate  quantities  of  powdered  titanium  metal  in  argon 
at  1000°C.  Dense  material  for  crucible  charges  was  then  obtained  by  arc-melting.  Melts  of 
both  compounds  could  be  contained  in  molybdenum  crucibles  with  little  contamination,  but  tanta¬ 
lum  crucibles  were  immediately  destroyed  by  melts  of  either  type.  Considerable  difficulty  was 
encountered  in  pulling  crystals  more  than  about  2  cm  long,  possibly  because  the  thermal  gradi¬ 
ent  in  the  solid  was  insufficient  to  stabilize  the  growing  interface. 

Use  of  an  induction  plasma  torch  to  grow  Ti^O^  by  the  Verneuil  method  is  also  being  studied. 
Large-grained  boules  have  been  obtained  in  the  initial  experiments. 

The  Seebeck  coefficient  a  of  a  Ti^O^  single  crystal  has  been  measured  between  7°  and 
550  °K.  The  results  of  the  measurements  are  shown  in  Fig.  Ill— 4 .  At  7°K,  a  =  —  2  p  volt /deg. 

With  increasing  temperature  a  becomes  positive,  increases  to  a  maximum  of  about  +270p.volt/ 
deg  at  210°K,  and  then  decreases  to  +1 2 \± volt/deg  at  550  °K.  A  somewhat  similar  positive  maxi- 

7 

mum  was  observed  by  Yahia  and  Frederikse,  who  measured  a  between  about  55°  and  1200°K. 

T.  B.  Reed  R.  E.  Fahey 
J.  M.  Honig  J.  T.  Roddy 


D.  CARRIER  CONCENTRATIONS  AND  SUPERCONDUCTIVITY  OF  In4SbTe3 

Hall  coefficient  measurements  have  been  made  at  77  °K  on  high-  and  low-pressure  samples 
of  In.SbTe-n  in  an  attempt  to  determine  whether  their  difference  in  superconducting  behavior  is 

**  s  8 

due  to  a  difference  in  carrier  concentrations.  As  previously  reported,  although  there  is  no  de¬ 
tectable  difference  in  structure  or  composition,  samples  prepared  at  high  pressure  exhibit 
superconductivity  at  temperatures  below  about  1.5  °K,  but  samples  prepared  at  atmospheric 
pressure  or  below  do  not  become  superconducting  down  to  1.3  °K,  the  lowest  temperature  at 
wiich  measurements  have  been  made. 

9 

The  Hall  coefficient  investigation  was  stimulated  by  a  recent  paper  by  Geller  and  Hull,  who 

proposed  an  ionic  model  to  account  for  the  metallic  behavior  and  superconductivity  of  high- 

+1  +3 

pressure  In^  ^Tefll).  According  to  this  model,  the  lattice  contains  In  and  In  ions  in  the 

ratio  required  to  balance  the  charge  on  the  Te  ^  ions.  Thus,  stoichiometric  InTe(II)  contains 

equal  numbers  of  In+1  and  In+^  ions,  while  indium-deficient  samples  contain  more  In+^  than  In+1 

ions.  The  high-pressure  phase  is  metallic  because  electrons  are  easily  transferred  from  one 

type  of  indium  ion  to  the  other.  The  concentration  of  free  electrons  is  equal  to  twice  the  con- 
11 

centration  of  In  ions  and  therefore  decreases  with  decreasing  indium  content.  According  to 
Geller  and  Hull,  this  decrease  in  electron  concentration  accounts  for  the  observed  decrease  in 
superconducting  transition  temperature  of  In^  xTe(II)  with  decreasing  indium  content. 
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Fig.  111-5.  Partial  pressures  of  Hg(g)  in  equilibrium  with  Hg-Te  condensed  phases, 
plotted  against  reciprocal  absolute  temperature.  Values  of  p^  along  three-phase 

lines  for  HgTe(c)  are  given  by  solid  curve,  those  for  Hg-Te  melts  by  dashed  lines. 


Fig.  1 1 1  —6.  Partial  pressures  of  Te^g)  in  equilibrium  with  Hg-Te  condensed  phases, 
plotted  against  reciprocal  absolute  temperature.  Values  of  Pje2  Ql°n9  three-phase 

lines  for  HgTe(c)  are  given  by  solid  curve,  those  for  Hg-Te  melts  by  dashed  lines. 
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On  the  basis  of  their  ionic  model,  Geller  and  Hull  predicted  and  prepared  a  number  of 
superconducting  phases  related  to  InTe(II),  including  In(Te^  ^As^  33)  and  In(Te0  ^As^  5>.  The 
analogy  between  these  two  materials  and  In^SbTe^,  which  Geller  and  Hull  did  not  consider,  sug¬ 
gested  that  the  ionic  model  if  valid  should  also  apply  to  In^SbTe^  and  that  differences  in  carrier 
concentration  might  explain  the  observed  differences  in  superconducting  behavior.  As  a  test  of 
th€;  model.  Hall  coefficient  (R^)  measurements  were  made  at  77  °K  on  three  low-pressure  and 
two  high-pressure  samples.  In  all  cases,  R^  is  negative,  indicating  conduction  by  electrons. 

The  carrier  concentrations  n  calculated  from  the  one-carrier  expression  n  =  — l/R^e  are  8.8, 

21  -3  22  -3  n 

9.0,  and  9.4  X  10  cm  for  the  low-pressure  samples  and  1.1  and  1.2  X  10  cm  for  the  high- 

pressure  samples.  The  latter  values  agree  rather  well  with  the  electron  concentration  of  1.3  X 
?2  -3 

10~  cm  calculated  for  In^SbTe^,  according  to  the  ionic  model.  Furthermore,  the  observed 

transition  temperatures  T  for  the  high-pressure  samples  are  1.4°  and  1.5  °K,  in  reasonable 

c  22  -3 

agreement  with  the  values  of  1.5°  to  1.7  °K  for  1.1  to  1.2  X  10  cm  ,  read  from  the  curve  of  n 
vs  T^  given  by  Geller  and  Hull  for  In^  xTe(II).  The  fact  that  low-pressure  samples  do  not  be¬ 
come  superconducting  down  to  1.3  °K  is  also  consistent  with  this  curve,  which  gives  T  =  1.1  °K 

21  -3  c 

for  n  =  9  X  10  cm  .  Thus,  the  model  proposed  by  Geller  and  Hull  accounts  for  several  fea¬ 
tures  of  the  data.  However,  since  the  model  predicts  that  n  is  determined  by  composition  only, 
it  cannot  explain  why  there  is  a  difference  in  concentration  between  the  high-  and  low-pressure 
samples,  which  apparently  have  the  same  composition.  In  order  to  make  a  further  test  of  the 
model,  it  is  planned  to  measure  the  Hall  coefficient  as  a  function  of  composition  for  samples  of 
In,  xTe(II)  and  for  high-pressure  InTe-InSb  alloys  containing  between  60  and  100  mol  %  InTe. 

M.  D.  Banus 
Lynne  B.  Farrell 
A. J.  Strauss 


E.  PARTIAL  PRESSURES  OF  Hg(g)  AND  Te2(g)  IN  Hg-Te  SYSTEM 


The  partial  pressures  of  Hg{g)  and  Te^fg)  in  equilibrium  with  Hg-saturated  HgTe(c),  Te- 
saturated  HgTe(c),  and  some  Hg-Te  melts  between  about  400°  and  700 °C  have  been  determined 
by  measuring  the  optical  density  of  the  coexistent  vapor  as  a  function  of  wavelength.  The  optical 
density  measurements  were  made  with  a  Cary  Model  14-H  double-beam  spectrophotometer,  as 
in  our  investigation  of  CdTe  (Ref.  10).  The  relationships  between  partial  pressure  and  optical 
density  were  obtained  from  the  results  of  calibration  experiments  on  pure  mercury  and  tellurium 
by  using  published  vapor  pressure  data  for  mercury11,12  and  tellurium.1^’ 14  Values  of 
were  calculated  from  the  measured  densities  at  a  number  of  wavelengths  between  2550  and 

o 


3 100 A.  At  these  wavelengths,  which  lie  on  the  long  wavelength  tail  of  the  mercury  absorption 
o  2 

line  at  2537  A,  the  optical  densities  were  found  to  be  proportional  to  (pH  )  ,  in  agreement  with 

less  extensive  data  reported  by  Kuhn  and  Freudenberg.  Values  of  PTe^  were  calculated  from 

the  measured  optical  densities  at  3650  and  4357  A,  which  are  proportional  to  PXe2'  ln  agreement 

with  mass  spectrographic  data,1^  the  absorption  spectra  give  no  evidence  for  the  existence  of 

HgTe(g)  molecules. 

The  measured  values  of  p^  and  PTe2  *n  equilibrium  with  condensed  phases  in  the  Hg-Te 
system  are  plotted  against  reciprocal  absolute  temperature  in  Figs.  Ill— 5  and  III— 6,  respectively. 
The  solid  lines  in  each  figure  are  three-phase  lines  which  give  the  limiting  values  of  partial 
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TABLE  lll-l 

LIQUIDUS  POINTS  IN  Hg-Te  SYSTEM 

Condensate 

Composition 

(atomic  percent  tellurium) 

Liquidus 

(°C) 

31.1 

638 

36.1 

645 

41.8 

659 

48.5 

670 

52.8 

666 

58.7* 

628 

*  Actual  tellurium  content  is  slightly  higher  be¬ 
cause  of  loss  of  mercury  during  preparation  of 
optical  cell. 

io3/t*k 


Fig.  Ill  —7.  Standard  Gibbs  free  energy  of  formation,  divided  by  absolute  temperature, 
for  reaction:  Hg  (g,  1  atm)  +  1/2  J&2  (g#  1  atm)  =  HgTe  (c,  equilibrium  vapor  pressure). 
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pressure  within  which  HgTe(c)  exists  as  a  stable  phase.  The  partial  pressures  of  both  elements 
are  seen  to  depend  strongly  on  the  composition  of  the  solid  phase.  Partial  pressures  in  equilib¬ 
rium  with  Hg-Te  melts  are  indicated  by  dashed  lines.  In  Fig.  Ill— 5,  the  upper  and  lower  branches 
of  the  three-phase  curve  give  values  of  p^  in  equilibrium  with  Hg-saturated  HgTe(c)  and  Te- 
saturated  HgTe(c),  respectively.  The  maximum  value  of  p^  in  equilibrium  with  HgTe(c)  is 
19  atm.  The  value  of  p^  at  the  maximum  melting  point  of  HgTe(c),  670°C,  is  about  12.5  atm. 

In  Fig.  III-6,  the  upper  and  lower  branches  of  the  three-phase  curve  give  values  of  Pxe^  eclu^" 

librium  with  Te-saturated  HgTe(c)  and  Hg-saturated  HgTe(c),  respectively.  The  maximum 

_  3 

value  of  PTez  equilibrium  with  HgTe(c)  is  3.1  X  10  atm. 

The  partial  pressure  data  for  Te-saturated  HgTe(c)  show  that  at  any  temperature  p^  is 
more  than  two  orders  of  magnitude  greater  than  PXe2*  Since  p^  has  its  minimum  value  and 
PTe2  has  *ts  maximum  value  for  Te-saturated  solid,  over  the  temperature  range  investigated 
Hg(g)  is  the  predominant  vapor  species  in  equilibrium  with  HgTe(c)  of  any  composition.  Thus, 
over  this  range,  HgTe(c)  has  no  congruently  subliming  composition,  and  continued  sublimation 
of  any  sample  will  cause  sufficient  loss  of  mercury  to  produce  Te-saturated  HgTe(c)  and  even¬ 
tual  formation  of  a  Te-rich  liquid  phase.  Extrapolation  of  the  experimental  results  shows  that 
this  is  also  true  down  to  at  least  300 °K. 

For  any  Hg-Te  condensate  containing  HgTe(c)  and  a  liquid  phase,  with  increasing  tempera¬ 
ture  the  values  of  p^  and  PTe2  along  the  three-phase  lines  of  Figs.  Ill- 5  and  III-6,  respec¬ 
tively.  At  the  liquidus  temperature,  where  the  condensate  becomes  entirely  liquid,  the  partial 
pressures  begin  to  rise  above  the  three-phase  lines.  If  the  pressure  change  is  sufficiently  ab¬ 
rupt,  the  liquidus  temperature  can  be  determined  from  the  intersection  of  the  three-phase  line 
wi*;h  the  partial  pressure  line  for  the  melt.  Liquidus  temperatures  obtained  in  this  manner  from 
the  data  of  Figs.III-5  and  III- 6  are  listed  in  Table  III—  1 .  The  condensate  compositions  given  have 
been  calculated  from  the  nominal  sample  compositions  by  correcting  for  the  preferential  loss  of 
mercury  to  the  vapor  phase. 

The  standard  Gibbs  free  energy  of  formation  of  HgTe(c),  in  equilibrium  with  the  vapor  phase, 

1  /2 

from  ideal  gaseous  Hg(g)  and  Te^(g),  each  at  1  atm,  is  given  by  AG°  [HgTe(c)]  =  RT  In  PpjgPTe2» 
if  it  is  assumed  that  the  vapor  phase  is  ideal  and  that  deviations  from  stoichiometry  are  negligible. 
Values  of  AG°/T  calculated  according  to  this  expression  from  the  experimental  data  for  p^^  and 
PTe2  are  Pl°ttec*  against  10^/T  in  Fig.  Ill- 7 .  The  estimated  errors  due  to  assuming  an  ideal 
vapor  phase  and  stoichiometric  solid  are  insignificant.  Most  of  the  values  of  AG°  were  obtained 
from  data  for  Te-saturated  HgTe(c);  those  obtained  from  data  for  Hg-saturated  HgTe(c)  are  in 
reasonably  good  agreement.  The  variation  of  AG°  between  778°  and  943  °K  is  well  represented 
by  the  linear  equation  AG°[HgTe(c)]  =—41.66  +  42.71  (10  Tkcal/mole.  Extrapolation  of  this 
equation  to  lower  temperatures  gives  values  of  AG°  which  agree  within  6  percent  with  those  ob¬ 
tained  by  Goldfinger  and  Jeunehomme^  from  effusion  measurements  between  435°  and  555 °K. 


R.  F.  Brebrick 
A.  J.  Strauss 
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F.  WET  CHEMICAL  ANALYSIS 
1.  HgTe 

A  method  has  been  developed  for  determining  mercury  and  tellurium  in  HgTe  with  an 

accuracy  of  about  one  part  per  thousand.  Mercury  is  determined  by  direct  automatic  titration 

1 7 

with  EDTA  to  a  potentiometric  endpoint,  using  the  Reilley  gold  electrode.  The  pH  is  main¬ 
tained  at  4.5  by  buffering  with  sodium  acetate-acetic  acid.  This  EDTA  titration  method  is  much 
faster  and  more  accurate  than  the  polarographic  method  previously  developed  for  mercury  in 

HgTe-CdTe  alloys.  Tellurium  is  determined  by  adding  excess  potassium  dichromate  and  back 

1 8 

titrating  with  ferrous  ammonium  sulfate  to  a  potentiometric  endpoint. 


2.  In-Sb-Te  System 

Methods  developed  for  determining  indium,  antimony,  and  tellurium  in  In-Sb-Te  samples 

19 

have  been  described  previously.  The  minimum  sample  size,  which  was  previously  0.5  gram, 
has  now  been  reduced  to  0.2  gram  by  the  use  of  more  dilute  titrants  (0.01  N)  in  determining 
indium  and  antimony.  No  reduction  in  titrant  concentration  could  be  made  in  the  tellurium 
determination,  since  inconsistent  results  were  obtained  with  weaker  solutions.  An  accuracy  of 
two  parts  per  thousand  in  determining  each  of  the  three  components  was  obtained  in  analyses  of 
synthetic  standards  weighing  0.2  gram.  j  ^  Qornwe^ 

E.  B.  Owens 
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IV.  BAND  STRUCTURE  AND  SPECTROSCOPY  OF  SOLIDS 


A.  INTERBAND  TRANSITIONS  IN  ANTIMONY* 

\ 

The  interband  transitions  in  antimony  which  have  been  reported  previously  have  now  been 
studied  in  more  detail.  A  summary  of  the  results  obtained  with  a  binary  face  in  the  photon  en¬ 
ergy  range  0.155  <ho;  <  0.305  ev  with  the  optical  electric  field  E  along  a  bisectrix  axis  is  given 
in  Fig.  IV- 1.  These  curves  can  be  experimentally  extrapolated  to  zero  field,  yielding  an  energy 
gap  of  0.143ev.  This  extrapolation  is  only  approximate,  since  the  lines  show  curvature,  im¬ 
plying  nonparabolic  energy  bands.  The  reduced  effective  mass  for  this  series  is  ~0.018  mQ, 
which  again  is  only  an  approximation  because  of  the  nonparabolic  character  of  the  bands.  Each 
oscillation  of  the  series  exhibits  structure  which  can  be  identified  with  two  almost  equal  effective 
masses,  or  with  a  spin  splitting.  If  the  former  interpretation  is  assumed,  then  the  two  reduced 
masses  are  ~0.0l6  mQ  and  ~0.019  mQ.  The  amplitude  of  these  oscillations  in  the  reflectivity  is 
very  small  for  E  along  the  trigonal  direction.  The  corresponding  interband  transitions  on  a 
bisectrix  face  have  also  been  studied.  The  reduced  effective  mass  here  is  approximately  equal 
;:o  that  for  the  binary  face  and  similar  polarization  effects  are  found. 

A  second  series  of  interband  transitions  is  observed  for  binary,  bisectrix,  and  trigonal  faces 
in  the  photon  energy  range  0.105  <  ho;  0.175ev.  Although  the  band  gap  associated  with  these  tran¬ 
sitions  (0.100  ±  0.003  ev)  is  smaller,  the  effective  masses  are  heavier  and  also  more  anisotropic. 

These  symmetry  considerations  suggest  that  the  two  sets  of  transitions  occur  at  different 
points  in  the  Brillouin  zone  (see  Fig.  IV-2)  with  the  larger  band  gap  occurring  at  point  T,  and  the 
smaller  band  gap  at  point  L.  A  more  detailed  account  of  this  work  has  been  submitted  for 

publication.  Mildred  S.  Dresselhaus 

J.  G.  Mavroides 


B.  ENERGY- MOMENTUM  DISPERSION  RELATIONS  FOR  GROUP  V  SEMIMETALS 

The  observation  of  two  series  of  interband  transitions  in  antimony  is  reported  in  Sec.  IV-A. 

2 

The  band  structure  calculation  of  the  Group  V  semimetals  indicates  that  interband  transitions 
are  likely  to  occur  about  the  points  L  and  T  in  the  Brillouin  zone.  We  present  here  the  dis¬ 
persion  relations  for  the  energy  bands  in  the  vicinity  of  these  two  points. 

Of  the  two,  the  point  L,  located  at  the  center  of  the  pseudohexagonal  face  (see  Fig.  IV-2), 
has  the  lower  symmetry,  which  is  designated  as  C^-  Using  IT  •  p  Brillouin-Wigner  perturbation 
theory^  and  the  symmetry  properties  of  this  point,  the  energy-momentum  relations  about  this 
point  can  be  derived.  The  notation  of  Ref.  2  is  used  to  label  the  bands  and  the  representations. 
The  wave  vector  k^  is  taken  along  the  binary  axis  and  transforms  as  the  representation  L^,  and 
:he  wave  vectors  k  and  k  ,  along  the  bisectrix  and  trigonal  directions,  respectively,  transform 

y  z  -*• 

as  the  representation  L^.  The  E  vs  k  relation  for  band  L.  is 


“This  experiment  was  carried  out  using  the  high  field  facilities  of  the  National  Magnet  Laboratory,  M.I.T. 
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Fig.  IV—  1 .  Summary  of  interband  transitions  for  a  binary  face.  Magnetic  field  is  perpendicular 
to  binary  face  and  optical  electric  field  is  along  a  bisectrix  axis.  T~  4°K. 


TRIGONAL  13-84-3705(01 


Fig.  IV-2.  Brillouin  zone  for  a  rhombohedral  crystal 
structure,  as  in  antimony.  Principal  symmetry  points 
and  directions  are  indicated. 
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ccuples  only  those  bands  given  in  Table  IV-1.  The  dispersion  relation  given  in  Eq.  (2)  holds 
equally  well  at  point  X,  the  center  of  the  rectangular  face. 

At  point  T,  located  at  the  center  of  the  real  hexagonal  face,  the  symmetry  follows  the  group 
D,^.  Here  there  are  four  one-dimensional  representations  (four  nondegenerate  bands)  and  two 
two-dimensional  representations  (two  sets  of  two  doubly  degenerate  bands).  The  wave  vector 
along  the  trigonal  direction  transforms  as  T^,,  while  the  wave  vectors  kx  and  k^  along  the  binary 
and  bisectrix  directions,  respectively,  transform  as  T^,.  The  E  vs  k  relations  for  the  four 
nondegenerate  bands  can  be  written  as 


I  z  |2 
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Table  IV-2  defines  the  momentum  matrix  elements  and  lists  the  selection  rules  for  the  terms 
in  these  summations. 


TABLE  IV-1 

SUMMARY  OF  SELECTION  RULES  AT  POINTS  L  AND  X 

i 

I 

L, 

L2 

L3 

L4 

L! 

0 

0 

X 

"1,3 

«Y'Z 

"1,4 

L2 

0 

0 

wy,z 

2,3 

X 

*2,4 

L3 

"3*. 

wy,z 

3,2 

0 

0 

L4 

X 

"4,2 

0 

0 
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TABLE  IV-2 

SUMMARY  OF  SELECTION  RULES  AT  POINTS  r  AND  T 
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The  energy  relations  for  the  and  T^,  bands  require  the  use  of  degenerate  k 
Wigner  perturbation  theory.  The  secular  equation 
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The  dispersion  relation  given  in  Eq.  (5)  holds  equally  well  at  point  T,  the  center  of  the  Brillouin 
zone. 

In  Eqs.  (1),  (3),  and  (5),  the  nonparabolic  character  of  the  bands  is  automatically  included 
in  the  k  dependence  of  the  energy  denominators.  Such  nonparabolic  effects  have  been  shown  to 
be  important  in  bismuth  (Ref.  4)  and  antimony  (Ref.  5).  The  effect  of  spin-orbit  interaction  can 
easily  be  included  in  this  formalism  and,  in  particular,  is  important  in  lifting  the  degeneracy 
of  the  T^  and  T^,  bands. 

These  dispersion  relations  are  being  applied  to  explain  the  magnetoreflection  experiments 

reported  in  Sec.  IV-A.  ^  ^ 

G.  Dresselhaus 

Mildred  S.  Dresselhaus 


C.  MAGNETOPLASMA  EFFECTS  IN  BISMUTH 

In  studying  the  magnetoplasma  effects  in  antimony  it  was  found  by  Dresselhaus  and 
Mavroides^  that  a  simple  three-tilted  ellipsoidal  model  failed  to  give  even  qualitative  agreement 
with  experiment.  Further  investigation  indicated  that  the  onset  of  direct  interband  transitions 
occurred  near  the  plasma  edges  suggesting  that  the  observed  incorrect  sign  and  amplitude  for 
the  magnetoplasma  effect  as  a  function  of  energy  for  certain  orientations  was  probably  due  to 
these  interband  effects.  In  order  to  establish  the  validity  of  the  free  multicarrier  model  calcu¬ 
lation,  experiments  are  now  being  carried  out  on  the  semimetal  bismuth  which  has  a  multi- 
el  Lipsoidal  Fermi  surface  somewhat  like  antimony  and,  in  addition,  the  further  simplification 
that  the  direct  interband  transitions  in  bismuth  occur  at  much  higher  photon  energies  than  the 
energies  corresponding  to  the  plasma  frequencies.  Preliminary  results  of  this  magnetoreflection 

7 

pr  eviously  reported  by  Boyle  and  Brailsford  are  not  detailed  enough  to  allow  a  comparison  with 
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Fig.  IV-3.  Variation  of  optical  reflectivity  with  photon  energy  for  H  =  3.  9  kgauss 
and  H  =7.  4  kgauss.  Light  is  incident  on  binary  face  and  magnetic  field  is  parallel 
to  a  bisectrix  axis. 


the  calculations.  Figure  IV-3  indicates  the  magnetoplasma  effects  measured  with  unpolarized 
light  incident  on  a  binary  face  with  the  magnetic  field  parallel  to  a  bisectrix  axis.  It  will  be  noted 
that  there  is  a  rather  complex  dependence  on  both  photon  energy  and  magnetic  field.  The  proper 
use  of  polarized  radiation  is  expected  to  allow  the  examination  of  each  plasma  edge  separately. 

D.  H.  Dickey 

D.  FARADAY  ROTATION  OF  INDIRECT  TRANSITION  IN  GERMANIUM* 

We  have  observed  on  transmission  the  oscillatory  Faraday  rotation  of  the  indirect  transition 
in  germanium.  These  preliminary  observations  were  made  at  magnetic  fields  of  62  kgauss  and 
at  helium  temperatures  on  a  heat  sunk  sample.  From  the  position  of  the  indirect  gap  the  temper¬ 
ature  was  determined  as  15°K.  The  sample  was  intrinsic,  with  a  carrier  concentration  n  < 

J  O  *5 

10  /cm  .  The  data  were  taken  with  the  magnetic  field  (and  the  direction  of  propagation)  perpen¬ 
dicular  to  a  (110)  face. 

Comparison  of  the  Faraday  rotation  with  the  corresponding  transmission  as  a  function  of 
wavelength  shows  that  there  are  oscillatory  curves  corresponding  to  both  the  exciton  absorption 
and  to  the  Landau  steps.  The  entire  indirect  transition  rotation  is  superposed  on  a  very  large 
dispersive  tail  which  probably  arises  from  the  r^,  —  T^,  direct  energy  gap  transition.  At  the 
temperatures  and  magnetic  fields  in  question  the  oscillatory  effects  are  of  the  order  of  2  percent 
of  the  total  rotation. 

Since  the  Faraday  rotation  is  a  dispersive  effect  and  the  direct  and  indirect  contributions  at 
a  given  wavelength  are  superposed,  one  should  be  able  to  make  a  direct  comparison  of  the  matrix 
elements  involved  in  the  respective  transitions.  Calculations  are  presently  being  carried  out 


toward  this  end. 


J.  Halpern 


*  This  experiment  was  carried  out  using  the  high  field  facilities  of  the  National  Magnet  Laboratory,  M.I.T. 
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E.  EXCITON  ABSORPTION  IN  GaSb’ 

g 

Exciton  lines  have  previously  been  observed  in  GaSb.  Some  of  these  have  been  attributed 
to  the  formation  of  complexes  consisting  of  an  ionized  acceptor  ion  and  an  exciton.  Similar  lines 
in  III-V  compounds  have  been  observed  in  emission  in  GaAs  (Sec.  IV-G)  and  InAs  (Ref.  9).  For 
such  a  complex  the  calculation  of  energy  levels  and  their  dependence  on  magnetic  field  appears 
to  be  rather  difficult.  The  study  of  the  absorption  lines  under  the  influence  of  magnetic  field  and 
under  the  influence  of  strain  may  provide  an  understanding  of  the  energy  levels  that  can  lead  to 
a  solution  of  the  theoretical  problem. 

Previous  worklu  has  shown  that  the  diamagnetic  effect  on  the  absorption  lines  (H  <  20kgauss) 
is  quadratic  and  is  the  same  as  that  on  the  free  exciton.  A  spin  splitting  is  observed  which  ap¬ 
pears  to  be  due  to  the  spin  of  the  electron  associated  with  the  complex.  Magnetic  field  measure¬ 
ments  on  the  y  line  (Ref.  8)  have  been  extended  to  ~90kgauss.  The  functional  behavior  of  the 
diamagnetic  shift  changes  beyond  20kgauss  and  eventually  becomes  linear.  The  critical  field 

ior  free  excitons  in  GaSb  is  ~20kgauss  and  the  observed  behavior  agrees  with  that  predicted  by 

1 1 

the  theory  of  Elliott  and  Loudon.  Analysis  of  the  observed  shift  yields  a  reduced  electron-hole 
effective  mass  of  0.051  mQ  which  may  be  compared  to  the  value  of  0.042  mQ  reported  previously. 
These  results  substantiate  the  earlier  indication  that  the  diamagnetic  effect  of  the  associated 
exciton-impurity  complex  is  the  same  as  that  of  the  free  exciton. 

The  observed  splitting  corresponds  to  |g|  =  10.3  which  is  slightly  higher  than  the  value 
of  | g  |  =  9.4  obtained  below  20kgauss  and  may  indicate  a  slight  increase  in  |g|  with  magnetic 
field. 

Changes  in  the  spectra  below  20kgauss  with  change  in  the  polarization  of  the  light  indicates 
structure  due  to  splitting  in  the  valence  band.  The  splittings  in  the  absorption  lines  in  fields  up 
1o  ~90kgauss  were  too  small  to  be  resolved.  ^  j  j0^nson 

G.  B.  Wright 


12 


F.  OPTICAL  ABSORPTION  BY  EXCITONS 

The  probability  per  unit  time  that,  with  a  radiation  density  p(v)  of  frequency  v ,  a  crystal 
will  make  an  electronic  transition  from  the  ground  state  4*^  to  an  excited  (exciton)  state  4'n#  K 
:.s  given  by 

W  =  — 1-J  |<*n'  K|  Yi  exp[iq  •  r.)  f  •  X  I*  >|2  PM  <5(E  -  E  +  hi')  (8) 

fiv  .  i  1  o  o 

where  q  and  £  are  the  wave  vector  and  polarization  vector  of  the  photon,  respectively,  and  j. 
:.s  the  current  operator  for  the  i  11  electron  at  the  position  r\.  In  order  to  evaluate  Eq.  (8)  we 
must  calculate  the  matrix  element  Z  exp[icf  •  F )  £  •  This  has  been  done  by 

13  1 

Elliott  in  the  limit  of  small  wave  vector  (q  =  0).  In  the  general  case  of  degenerate  conduction 
and  valence  bands  the  result  in  this  limit  is 


"This  experiment  was  carried  out  using  the  high  field  facilities  of  the  National  Magnet  Laboratory,  M.I.T. 


41 


Section  IV 


s  t 


<*n|  Z  r  •  ^  l*0>  =  V1/2  Yj  E  *cnv.(r)  I  +  [_iVr  ^c.v.^*]  I  Vk 

i  i=l  j=l  i  r=0  i  ^  r=0 

+  f  f(-iVa)  (-iVP)  (i 

2  1'  r  r  Tc.v. 

1  i  J 


V.a  vf 
k  k 


+  .  .  . 


<*c  kl  { 

i 


pv.k/ 

J 


k=ku 


(9) 


where  i  and  j  run  over  the  s  degenerate  conduction  bands  and  t  degenerate  valence  bands, 
respectively,  $c.v,(r)  is  a  solution  of  a  set  of  s  times  t  simultaneous  differential  equations 
(exciton  effective  mass  equations;*"*  r  is  the  relative  coordinate  of  the  electron  and  hole  and 
Eq.  (9)  is  evaluated  at  r  =  0;  a  and  /3  run  over  the  directions  x,  y,  and  z  and  the  operators 

take  the  derivative  of  the  band-to-band  matrix  element  £  *  j  l^vjk^w^  resPeci  to  a 

variation  of  ka  evaluated  at  k  =  k°  the  location  in  k-space  of  the  direct  band-to-band  transition. 

For  finite  wave  vector  q  and  in  the  presence  of  external  magnetic  fields  we  have  found  the 
generalization  of  Eq.  (9)  to  be, 
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where  the  operator  D  is  given  by 
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where  V  takes  the  derivative  of  if  **' *v  with  respect  to  K.  t/in’^  is  an  implicit  function  of  K  since 

_  ciVj  ciVj 

the  Hamiltonian  of  which  ip  is  an  eigenfunction  depends  on  K  (Ref.  14).  In  addition,  the  only 

CiVj 

difference  between  Eq.  (10)  and  Eq.  (9)  is  that  whereas  in  Eq.  (9)  the  conduction  and  valence  band 
functions  were  both  evaluated  at  k  =  k°  in  Eq.  (10),  the  conduction  and  valence  band  functions  are 
evaluated  at  k  =  k°  +  l/2q  and  k  =  k°  —  l/2q,  respectively,  in  terms  of  an  expansion  about  k°. 
That  is,  we  write 


*  O  ,  =  exp[i(k  +  2  q)  •  r  ] 
c.,k°+q/2  Z  uc.,k°+q/2 


and 
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sjid  similarly  for  ip  .  Therefore,  the  matrix  element  ^1  2  expfiqT .  F.l  £  •  7.  1^  ^ 

^  vj,k°-q/2  ^  '  i  l  1  '  ° 

can  be  easily  obtained  as  an  expansion  in  powers  of  q  from  Eq.  (10). 

J.  O.  Dimmock 

G.  PHOTOLUMINESCENCE  OF  GaAs  IN  A  HIGH  MAGNETIC  FIELD* 

The  photoluminescence  of  GaAs  has  been  studied.  The  lines  observed  in  n-type  material 

15 

included  those  reported  at  zero  field  by  Nathan  and  Burns.  The  exciton  lines  and  several  tran¬ 
sitions  proceeding  through  holes  trapped  at  acceptors  all  exhibited  the  same  variation  of  photon 
energy  with  magnetic  field.  The  energy  shift,  which  was  quadratic  at  low  fields  and  linear  at 
high  fields,  could  be  described  by  a  simple  effective  mass  hydrogenic  model  with  an  effective 
mass  for  the  electron  of  about  0.06  mQ,  and  a  binding  energy  of  approximately  5  mev.  The 
binding  energy  required  for  the  fit  changed  very  little  with  the  zero-field  photon  energy  of  the 
emission  lines,  which  in  some  cases  was  as  much  as  200  mev  below  the  energy  gap  of  GaAs. 
Evidently,  the  binding  energy  of  the  trapped  hole  had  little  effect  on  the  magnetic  behavior  of  the 
3*ecombination  radiation.  Attempts  to  resolve  a  spin  splitting  due  to  holes  were  unsuccessful  for 
any  of  the  transitions.  The  g-factors  for  holes  on  zinc  and  cadmium  impurities  reported  by 
Title ^  for  strained  samples  would  have  given  rise  to  splittings  well  within  the  resolution  capa¬ 
bilities  of  the  instrument.  In  fact,  a  fourfold  splitting  of  the  exciton  line  in  CdTe  has  been 
observed. 

In  the  n-type  material  studied,  the  intensity  of  the  lines  was  practically  unaffected  by  the 
magnetic  field.  In  p-type  material,  however,  an  additional  line  was  observed,  and  as  the  mag¬ 
netic  field  increased,  this  line  shrank,  while  the  other  lines  grew,  sometimes  as  much  as  ten¬ 
fold  in  intensity.  Evidently,  the  field  was  changing  the  relative  electron  capture  cross  sections 
for  the  various  centers  involved.  The  magnitude  of  the  effect  was  not  reproducible  with  time. 

G.  B.  Wright 
F.  L.  Galeener^ 

H.  PIEZOELECTRIC  POLARONS  IN  STRONG  MAGNETIC  AND  COULOMB  FIELDS 

17  18 

Recently  reported  measurements  '  of  the  cyclotron  resonance  frequency  of  electrons  in 

the  conduction  band  of  CdS  seem  to  indicate  that  the  Landau  levels  are  approximately  20  percent 

more  widely  separated  than  is  to  be  expected  on  the  basis  of  the  electronic  effective  mass.  A 

theory  which  purports  to  explain  this  anomaly  on  the  basis  of  piezoelectric  electron-phonon 

19 

interaction  has  been  proposed  by  Hopfield  and  Mahan,  but  this  theory  is  subject  to  grave  crit¬ 
icism  on  theoretical  grounds. 

Like  Hopfield  and  Mahan,  we  have  used  a  Frohlich-type  Hamiltonian  with  an  isotropic 
electron-phonon  interaction  term  to  serve  as  a  model  for  the  piezoelectric  electron-phonon 
coupling.  However,  we  have  studied  the  effects  of  such  coupling  on  the  electron  energy  at  zero 
temperature  (where  the  lowest  order  perturbation  theory  can  be  rigorously  evaluated)  in  the 
presence  of  (1)  a  uniform  external  magnetic  field  H,  and  (2)  a  coulomb  potential  binding  the 
electron  to  a  fixed  center. 


*This  experiment  was  carried  out  using  the  high  field  facilities  of  the  National  Magnet  Laboratory,  M.I.T. 
t  National  Magnet  Laboratory,  M.I.T. 
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In  CdS,  the  dimensionless  piezoelectric  coupling  constant  a ,  defined  so  as  to  play  the  same 

role  as  Frohlich's  a  plays  in  the  usual  polaron  theory,^  is  found  to  be  approximately  0.2  (using 
21 

Hutson's  data). 

Therefore,  we  have  treated  the  electron-phonon  interaction  in  lowest-order  perturbation 

theory  in  the  presence  of  a  strong  magnetic  field  in  order  to  determine  the  perturbation  of  the 

energy  spacing  between  the  n  =  1  and  n  =  0  Landau  levels  at  zero  temperature.  In  the  limit  of 

1  /4 

infinite  field,  the  spacing  is  found  to  increase  by  an  amount  proportional  to  H  '  .  However,  for 
magnetic  fields  readily  attainable  in  the  laboratory,  it  is  necessary  to  evaluate  the  perturbation 
expression  numerically.  Preliminary  results  indicate  that  for  field  strengths  of  the  order  of 
several  kilogauss  the  piezoelectric  interaction  in  this  model  is  much  too  weak  to  account  for  the 
20-percent  effect  observed. 

If  we  denote  the  discrete  index  which  labels  the  energy  of  an  electron  in  a  magnetic  field  by 
n,  then  we  find  that  the  "golden  rule"  expression  for  the  decay  probability  of  electron  levels, 
which  can  emit  a  phonon  without  change  of  n  and  still  conserve  energy  in  the  process,  diverges. 
Closer  analysis  shows  that  lowest-order  perturbation  theory  implies  that  the  decay  is  nonexpo¬ 
nential  in  time.  The  consequences  for  the  line  shift  at  temperatures  above  zero  need  further 
investigation. 

The  effect  of  the  piezoelectric  electron-phonon  interaction  on  the  ionization  energy  of  an 

electron  bound  in  a  coulomb  field  to  a  point  center  of  charge  is  to  increase  the  ionization  energy 

by  a  percentage  which  is  independent  of  the  strength  of  the  coulomb  field  for  strong  fields.  In 

CdS  this  should  be  approximately  a  2-percent  effect. 

Interesting  results  for  upper  and  lower  limits  of  the  ground  state  energy  of  a  polaron  bound 

in  a  coulomb  field  have  been  obtained.  n/I  T 

L).  M.  Larsen 


I.  CALCULATION  OF  RELATIVISTIC  ENERGY  BANDS  IN  CRYSTALS 


22 


In  connection  with  a  continuing  interest  in  the  band  structure  of  rare  earth  metals  we  have 

obtained  a  formal  solution  of  the  relativistic  Schrodinger  equation  for  the  periodic  "muffin  tin" 

23 

potential  in  the  framework  of  the  augmented  plane-wave  (APW)  method.  The  problem  of  eval¬ 
uating  the  relativistic  solution  for  specific  materials  does  not  appear  to  be  much  more  difficult 
than  that  of  evaluating  the  nonrelativistic  solution.  The  principal  inaccuracies  in  the  relativistic 
calculation  are  still,  as  in  the  nonrelativistic  calculation,  due  to  inaccuracies  in  the  assumed 
effective  potential.  Disregarding  these  errors,  it  appears  possible  to  obtain  solutions  of  the 
relativistic  periodic  potential  problem  to  an  accuracy  of  about  ±0.001  with  present  day  com¬ 
puters.  This  would  result  in  the  calculation,  directly  using  an  APW  method,  of  spin-orbit 
splittings  and  other  relativistic  effects  in  solids  directly.  The  method  appears  quite  feasible 
and  attention  is  now  being  given  to  the  problems  of  programming  the  calculation  to  obtain  rela¬ 
tivistic  energy  bands  in  some  specific  materials.  j  q  p)^mmocj< 


J.  OPTICAL  STUDIES  ON  SULFUR- DOPED  SILICON 

24 

We  have  been  studying  the  optical  properties  of  sulfur-doped  silicon.  By  suitably  choosing 
the  starting  material  and  thermal  treatment,  one  can  obtain  samples  with  sulfur  impurities  ex¬ 
isting  as  neutral  single  sulfur  atoms,  singly  ionized  single  sulfur  atoms,  or  as  singly  ionized 
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sulfur  pairs.  The  absorption  lines  due  to  sulfur  in  silicon  are  shown  in  Fig.  IV-4,  where  a  bracket 

indicates  a  group  of  lines  due  to  a  given  type  of  sulfur  center.  When  these  absorption  lines  were 

studied  under  stress  it  was  found  that  all  the  absorptions  due  to  transitions  between  a  ground 

state  and  the  excited  p-like  states  had  identical  shifts.  Figure  IV-5  shows  the  splitting  and  the 

shift  of  energies  of  the  transitions  to  p-like  states  under  an  applied  stress  in  the  [001]  direction 

for  each  type  of  sulfur  center.  Measurements  were  also  made  using  polarized  light.  Figure  IV-6 

shows  the  absorption  spectra  of  the  neutral  single  sulfur  atoms  at  zero  stress,  and  at  a  stress  of 
8  Z 

2.24  x  10  dynes/cm  for  light  polarized  both  parallel  and  perpendicular  to  the  applied  stress. 

For  a  uniform  longitudinal  stress  T  =  Tmm,  where  m  is  a  unit  vector  along  the  stress 

axis  with  components  m  ,  m  ,  and  m  referred  to  the  cube  axes  of  the  crystal,  the  minimum 
th  x  y  z 

of  the  j  valley  of  the  conduction  band  in  silicon  shifts  by 

6E(j)  =  E 

I'll  1£  £11 
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1(s11  +  2s12)T  +  E2(sH  Sl2)  (m|j|  3)T 


(14) 


where  j  =  x,  x,  y,  y,  z,  z;  E7  is  Brook’s  notation  for  the  pure  shear  deformation  potential 

z  26  *-* 

coefficient  of  a  valley,  which  is  denoted  by  Herring  and  Vogt  by  and  s^  and  s ^  are  the 
standard  elastic  compliance  coefficients.  Equation  (14)  assumes  that  the  shift  in  energy  of  a  valley 
minimum  depends  linearly  on  the  strain  (or  the  stress).  In  the  effective  mass  approximation, 
strain  will  shift  the  donor  levels  in  this  manner  when  intervalley  coupling  is  negligible.  When 
T  is  parallel  to  the  [001]  direction,  the  equation  for  the  splitting  of  a  p-state  reduces  to  the  form 

6E(z,  z)  _  6E(x>x,y,y) 


E2  = 


<s  1 1  s12)  T 
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Fig.  IV-4.  Absorptions  in  sulfur  doped  silicon  for  samples  with  various  starting  materials 
and  thermal  treatments. 
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(S~S)+  IMPURITY  1  3-H-36TI(u] 


S+  IMPURITY 


STRESS  (x  10® dynes/cm2) 

Fig.  IV-5.  Effects  of  stress  on  transitions  to  p-like  states  for  three  types  of  sulfur 
impurity  centers. 


Fig.  IV-6.  Effect  of  stress  on  S°  impurity  spectrum  in  silicon. 
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Using  this  equation,  and  substituting  the  measured  value  for  the  splitting  of  the  p-like  states 

27  _  j  2 

urder  stress,  and  Fine's  value  of  (s^  —  s^)  =  1.016  dynes/cm  ,  we  obtain  the  value  E 2  = 

7.9  ±  0.2  ev. 

W.  E.  Krag  H.  J.  Zeiger 
W.H.  Kleiner  S.  Fischler 


K  UNUSUALLY  LARGE  ELECTRON  HALL  MOBILITY,  MAGNETORESISTANCE 
AND  MAGNETO-HALL  EFFECTS  IN  HgTe 

Hall  coefficients  R  and  magnetoresistance  Ap/pQ  have  been  measured  at  4.2°  and  78 °K  on 
a  variety  of  HgTe  samples  in  magnetic  fields  H  ranging  from  0.4  to  105kgauss.  The  dependence 
of  R  on  H  for  a  selected  set  of  samples  is  shown  in  Fig.  IV-7,  and  the  key  data  pertaining  to 
these  and  other  samples  are  summarized  in  Table  IV-3  in  order  of  decreasing  estimated  ratios 
of  electron  n  to  hole  p  concentrations.  Sample  A  is  typical  of  an  extrinsic  n-type  material 
below  78 °K.  Corresponding  to  a  large  n/p  ratio,  R  is  negative  and  independent  of  H  and 
Ap/p  ~  1.  As  one  proceeds  to  smaller  n/p  ratios  in  samples  B  through  E,  the  low  field  Hall 
coefficients  R(o)  pass  through  a  deep  minimum  of  about  —3000  cm  /coulomb,  and  Ap/po  »  1. 

Trie  large  numerical  values  of  R  at  4.2 °K  in  low  magnetic  field  are  consistent  with  the  overlap 
band  model  of  HgTe  (Ref.  28),  since  the  electron  to  hole  mobility  ratios  are  large  and  n/p  <  1. 
ALso,  for  the  range  of  n/p  values  encountered  in  samples  D  to  F,  R(H)  at  4.2 °K  changes  from 
negative  to  positive  with  increasing  H.  This  behavior  is  expected  for  a  material  with  an 


Fig.  IV-7.  Experimental  curves  of  Hall  coefficient  R  vs  magnetic  field  strength  for  samples 
B,  D,  F,  and  G  of  HgTe.  Solid  curves  designate  negative  Hall  coefficients,  whereas  dashed 
curves  denote  positive  R.  (High  field  measurements  were  made  at  the  National  Magnet  Labo¬ 
ratory,  M.l.T.) 
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overlapped  band  model.  Sample  G  is  typical  of  extrinsic,  p-type  material;  here  n/p  —  0,  R  >  0 
for  all  H,  and  A p/pQ  <  1. 

It  is  of  interest  that  samples  for  which  R  >  4  cm3/coulomb  have  not  been  encountered.  This 
is  in  agreement  with  theoretical  predictions  which  place  an  upper  limit  on  positive  R  for  semi¬ 
metals  with  large  electron  to  hole  mobility  ratios. 

The  results  cited  can  be  qualitatively  understood  as  follows:  In  sample  A  the  Fermi  level 
Ep  is  at  least  4kT  above  the  valence  band  (VB)  edge;  hence,  only  electrons  contribute  to  trans¬ 
port  properties.  Sample  B  is  nearly  intrinsic;  here,  lies  below  the  top  of  the  VB,  so  that 
both  holes  and  electrons  participate  in  conduction,  rendering  A p/pq  ~  200  at  105kgauss.  In 
passing  from  sample  C  to  G,  Ep  systematically  diminishes  until  it  is  more  than  4kT  below  the 
bottom  of  the  conduction  band  edge.  At  4.2°K,  only  holes  conduct  in  sample  G;  however,  the 
thermal  spread  in  electron  energies  at  78 °K  is  sufficient  to  enable  high  mobility  electrons  to 
participate,  thus  making  R(o)  <  0  at  this  temperature.  Positive  R  were  encountered  both  at 

4.2°  and  78 °K  for  samples  more  heavily  doped  by  p-type  impurities  than  sample  G. 

2 

Table  IV-3  shows  that  electron  Hall  mobilities  in  excess  of  350,000  cm  /volt-sec  have  been 
observed  at  4.2°K.  These  values  are  higher  than  any  mobilities  reported  for  HgTe  to  date  and 
are  in  excess  of  the  mobilities  cited  for  any  zinc  blende  material,  including  InSb,  at  this  tem¬ 
perature.  These  high  electron  mobilities  are  not  due  to  the  large  static  dielectric  constant  as 
in  PbTe  but  are  associated  with  the  low  electron  effective  mass  and  the  fairly  high  state  of  purity 
of  the  material  under  study. 

Attempts  to  fit  these  data  quantitatively  to  a  mixed  conduction  model  are  in  progress. 

T.  C.  Harman 
J.  M.  Honig 
A.  E.  Paladino 


L.  THERMAL  ENERGY  GAP  OF  HgTe 

The  thermal  energy  gap  of  HgTe  has  been  determined  from  data  on  the  Hall  coefficient  R  by 

the  same  method  previously  used  to  find  the  thermal  energy  gap  of  HgSe  (Ref.  29).  The  data  were 

obtained  in  experiments  on  n-type  samples  with  net  donor  concentrations  between  2  x  10  ^  and 
18  —  3 

5X10  cm  as  determined  from  Hall  coefficient  measurements  at  4.2°K.  (By  definition,  a 
material  is  n-type  if  n  >  p  and  p-type  if  p  >  n.)  These  samples,  which  were  cut  from  large¬ 
grained  ingots  grown  from  iodine -doped  melts  by  the  Bridgman  technique,  were  mercury  satu¬ 
rated  by  annealing  at  350°C  to  reduce  the  concentration  of  acceptor  defects  resulting  from  stoi¬ 
chiometric  deviations.  Identification  of  samples  as  n-type  was  based  on  values  of  resistivity 
and  Hall  coefficient  measured  at  300°,  77°,  and  4°K  for  various  values  of  magnetic  field  strength. 
This  detailed  identification  was  necessary  because  some  p-type  samples  of  HgTe  have  negative 
Hall  coefficients  even  at  4°K,  since  the  ratio  of  electron  mobility  to  hole  mobility  is  at 
least  100  in  HgTe  (Ref.  30).  Because  of  this  high  mobility  ratio,  the  electron  concentration  in 
an  intrinsic  or  n-type  sample  at  any  temperature  is  given  by  the  usual  expression  for  one-carrier 
conduction:  n  =  —  l/Rec.  For  all  n-type  samples,  the  electron  concentration  either  remained 
essentially  constant  or  first  remained  constant  and  then  increased  monotonically  with  increasing 
temperature  from  4°  to  300°K.  The  increase  in  concentration  is  due  to  thermal  excitation  of 
electrons  from  a  valence  band  rather  than  from  a  partially  ionized  donor  level. 
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TABLE  IV-3 

SUMMARY  OF  PERTINENT  DATA  FOR  SAMPLES  A  TO  G 


Sample 

No. 

Sample 

Type 

R4(°)* 

P4(0)f 

r78(°)* 

‘W0*’ 

MH78 

A 

N 

-40 

5.3X10"4 

7.6X104 

«1 

-40 

5.6  X10"4 

7.1  XlO4 

B 

P 

-3000 

8.4X10'3 

3.6X105 

200 

-98 

1.2X  10"3 

8.  2  XlO4 

C 

P 

-1500 

1.5X10"2 

1  xio5 

- 

-47 

1.1  xio-3 

4.  3  XIO4 

D 

P 

-2800 

3.3X10-2 

3. 6  X  104 

27 

-99 

2.1  XIO"3 

4. 7  XIO4 

E 

P 

-960 

1.2X  10"2 

8  XIO4 

— 

-66 

7. 5  X  1 0-3 

8. 8  XIO3 

F 

P 

-5 

1.4X  10-2 

2.9  XlO2 

<1 

-15 

1 . 3  X lO-2 

1.2  XIO3 

G 

P 

+  1.1 

2.1  X10-2 

52 

<1 

-5 

1.6  XlO"2 

3.1  XIO2 

*Weak  field  Hall  coefficient  at  4°  and  78°K  in  cm  /coulomb. 

t  Electrical  resistivity  for  zero  magnetic  field  at  4°  and  78°K  in  ohm-cm. 

2 

$  Hall  mobility  at  4°  and  78° K  in  cm  /volt-sec. 

§  Magnetoresistance  for  a  magnetic  field  of  105  kgauss  at  4°K. 
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Fig.  IV-8.  Experimental  and  theoretical  values  of  electron  concentration  at  300°K(ngoo) 
vs  electron  concentration  at  77° K(n^)  for  HgTe.  Theoretical  values  were  calculated  using 

m^^=0. 16#  E~=0.14ev,  P  =6.5X10  ^  ev-cm,  and  the  values  of  E  given  in  the 

p  kj  t 

figure. 


The  initial  theoretical  plot  of  n~nn  vs  n77  did  not  agree  satisfactorily  with  experiment,  since 
(N)  J 

the  values  of  E,  and  m  were  only  preliminary.  In  order  to  obtain  better  agreement,  the  pro- 
1  P  (n) 

cedure  described  above  was  repeated  with  different  combinations  of  E,  and  m  values.  The 

(N)  t  p 

best  fit  to  the  data  was  obtained  for  m  =  0.16  m  .  The  degree  of  agreement  is  shown  in 

P  °  (fsj) 

Fig.  IV-8,  which  gives  the  experimental  results  and  theoretical  curves  for  m^  =  0.l6mQ  and 
three  values  of  E^.  Because  of  the  sensitivity  of  the  calculated  results  to  the  value  of  E^,  it  can 


be  concluded  that  E+  =  —  0.02  ±  0.01  ev.  For  HgSe,  values  of  E  =  —0.07  ev  and  m 

1  29  1 

were  obtained  previously. 


(N) 


=  0.17m 


T.  C.  Harman 
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V.  MAGNETISM  AND  RESONANCE 


A.  THEORETICAL 

1.  Determination  of  Ordering  Temperatures  from  High- Temperature  Expansions 

Kramers  and  Opechowski  first  suggested  that  one  could  estimate  the  Curie  temperature 

from  the  high-temperature  expansion  for  the  zero-field  susceptibility  per  particle  x(T)  by  using 

the  definition  that  x(T)  00  as  T  -*•  T  .  The  susceptibility  can  accordingly  be  expressed  as  a 

-1  c 

power  series  in  fi  =  (kT)  ,  or  more  conveniently  as 


X(T)  =  l  arrr 
r=0 


(1) 


T^/T,  T^  being  the  ordering  temperature  estimated  by  the  molecular-field  approxi- 


where  r 

mation.  Given  a  particular  model,  the  first  few  terms  of  this  series  can  be  evaluated  explicitly, 
and  an  estimation  of  its  radius  of  convergence  can  be  attempted.  If  the  ratio  test  is  used  to  ob¬ 
tain  this  estimate,  one  finds  that  T  =  pT**,  where  p  =  lim  la  ,  ./a  I  .  If  the  ratios  of  suc- 

c  M  1  n+1  n1 

n-*oo 

cessive  pairs  of  terms  are  sufficiently  well  behaved,  this  estimate  can  be  made  with  reasonable 
accuracy;  otherwise  it  cannot. 

We  have  carried  out  the  necessary  calculations  for  the  case  of  a  normal,  cubic  spinel  pos¬ 
sessing  both  A-B  and  B-B  nearest -neighbor  exchange  interactions.  We  used  the  quantum- 
rcechanical  Heisenberg  Hamiltonian  and  obtained  x(T)  to  fourth  order  in  r.  The  computed  ratios 
are  listed  in  Table  V-l  for  various  values  of  the  exchange  parameter  u  = 

This  method  for  estimating  T  is  essentially  that  which  Wojtowicz2*  has  applied  to  the  case  of 
spinels  with  only  A-B  interactions,  and  indeed  the  successive  ratios  shown  in  Table  V-l  for 
u  =  0  appear  to  be  reasonably  well  behaved.  However,  they  become  unmanageable  when  J-g-g  be¬ 
comes  appreciable,  whether  it  is  ferromagnetic  (u  <  0)  or  antiferromagnetic  (u  >  0).  It  is  gen¬ 
erally  believed  that  the  molecular  field  approximation  overestimates  the  ordering  temperature, 
usually  by  an  appreciable  amount.  Thus,  the  results  listed  in  Table  V-l  for  positive  u  preclude 

any  possibility  of  obtaining  even  a  crude  approximation  to  T  by  this  method. 

3  c 

In  a  recent  paper,  we  introduced  the  function  defined  by 


<p(c,  T)  =  N-1  2  c?c.  <S.  •  S.>T/(S.S.) 

ij 


(2) 


_ 2  th 

wnere  S ?  =  +  1),  Si  being  the  spin  quantum  number  for  the  1  site,  and  c  denotes  a  set  of 


c^,  which  are  constrained  by 

Z  I  I  2  =  N 


(3) 


Given  the  high-temperature  expansion  for  the  normalized  spin  correlation  function  ^ 


(S.S.),  one  can  write 

1  j" 
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TABLE  V-l 

RATIOS  OF  SUCCESSIVE  COEFFICIENTS 

IN  POWER  SERIES  REPRESENTATION 

MAGNETIC  SUSCEPTIBILITY  FOR  S.  =  S.  = 

A  B 

=  3/2 

u 

o 

< 

o 

CN 

O 

°3/o2 

°4/a3 

2.00 

-2.743 

-2.485 

-2.349 

-2.287 

1 .60 

-2.279 

-2.053 

-1.959 

-1.907 

1.20 

-1.870 

-1.683 

-1.614 

-1.570 

0.80 

-1.513 

-1.372 

-1.308 

-1.277 

0.40 

-1.205 

-1.118 

-1.031 

-1.031 

0.00 

-0.943 

-0.919 

-0.772 

-0.849 

-0.40 

-0.721 

-0.775 

-0.513 

-0.791 

-0.80 

-0.535 

-0.691 

-0.239 

-1.219 

-1.20 

-0.379 

-0.681 

0.056 

5.261 

-1.60 

-0.250 

-0.791 

0.346 

1.103 

-2.00 

-0.142 

-1.191 

0.581 

0.856 

TABLE  V-2 

RATIOS  OF  SUCCESSIVE  COEFFICIENTS  IN  POWER  SERIES 
REPRESENTATION  OF  CORRELATION  OBTAINED 

WITH  MOLECULAR-FIELD  SPIN  CONFIGURATIONS  (S  =S  =3/2) 

u 

a 

> 

o 

°2/al 

CN 

CO 

o 

a4/a3 

2.00 

1.000 

0.493 

0.541 

0.046 

1.60 

1.000 

0.627 

0.597 

0.443 

1.20 

1.000 

0.732 

0.668 

0.657 

0.80 

1.000 

0.810 

0.732 

0.770 

0.40 

1.000 

0.864 

0.784 

0.829 

0.00 

1.000 

0.899 

0.821 

0.858 

-0.40 

1.000 

0.917 

0.845 

0.869 

-0.80 

1.000 

0.925 

0.857 

0.869 

-1.20 

1.000 

0.924 

0.860 

0.863 

-1.60 

1.000 

0.918 

0.857 

0.853 

-2.00 

1.000 

0.910 

0.850 

0.840 

54 


Section  V 


oo 


<p{c,  T)  =  Yi  arT 


r 


(4) 


r=0 


in  analogy  with  Eq.  (1).  Each  c  (that  is,  set  of  c.)  represents  a  spin  configuration,  and  the  a 
of  Eq.  (4)  depend  upon  this  c.  In  effect,  the  function  <p( c,  T)  serves  to  project  the  spin  correla¬ 
tions  onto  the  configuration  c.  For  example,  if  c  is  chosen  so  as  to  correspond  to  a  ferro¬ 
magnetic  spin  arrangement,  then  Eq.  (4)  becomes  equivalent  to  Eq.  (1). 

However,  we  are  concerned  with  finding  the  correction  to  the  molecular-field  ordering 
temperature  for  a  class  of  ferrimagnetic  materials.  In  such  a  context,  it  would  appear 
more  reasonable  to  choose  the  c  so  as  to  represent  the  ferrimagnetic  spin  configuration  com¬ 
puted  in  the  molecular-field  approximation  for  T  s  T^.  We  have  carried  out  this  calculation 
up  to  order  r  =  4,  and  the  ratios  of  the  resulting  a^  are  listed  in  Table  V-2  for  the  same  set  of 
u-values  as  before.  We  note  first  that,  by  virtue  of  the  above  construction,  the  initial  ratio 
a^/a q  must  always  be  unity.  Then  we  note  the  general  improvement  of  the  behavior  of  suc¬ 
cessive  ratios  over  that  presented  in  Table  V-l.  There  is  significant  improvement  for  u  =  0; 
the  improvement  for  nonzero  u  is  pronounced.  These  new  results  are  quite  satisfactory  for 
negative  u,  where  the  A-B  and  B-B  interactions  assist  each  other.  Nevertheless,  they  deteri¬ 
orate  for  sufficiently  large,  positive  u,  where  the  interactions  compete  strongly. 

In  the  previous  Solid  State  Research  Report,  we  were  primarily  concerned  with  the  problem 
of  estimating  the  type  of  magnetic  ordering  to  be  expected  at  Tc  from  the  behavior  of  c,  where 
c(T)  was  defined  as  that  c  which  maximizes  <p{ c,  T)  (Ref.  3).  In  other  words,  the  configuration 
c(T)  possesses  the  largest  possible  projection  of  the  spin  correlations.  The  function  <p(c,  T)  can 
then  be  expressed  as 


oo 


(5) 


r=0 


which  differs  from  the  power  series  of  Eq.  (4)  in  that  the  implicit  dependence  of  the  original  a^ 
upon  T  through  their  dependence  upon  the  varying  configuration  c  (T)  has  been  explicitly  taken 
into  account.  Since  <p{ c,  T)  is  the  maximum  </?( c,  T),  it  must  diverge  the  most  strongly  at  Tc- 
Hence  the  estimation  of  a  radius  of  convergence  should  follow  most  readily  from  the  series 
given  in  Eq.  (5) . 

The  coefficients  br  have  been  computed  up  to  order  r  =  4,  and  the  resulting  ratios  are  listed 
in  Table  V-3.  The  first  two  columns  are  identical  with  those  of  Table  V-2  because  the  molecular- 
field  solution  constitutes  the  first  approximation  to  c  (T),  namely  its  temperature -independent 
tcirm  c  («).  However,  the  improvement  in  the  behavior  of  the  higher  ratios  is  appreciable  even 
for  u  =  0,  and  becomes  large  for  very  large  u.  For  example,  the  modifications  in  the  sequence 
of  ratios  for  u  =  1.6  are  particularly  striking,  as  one  compares  the  different  tabulations.  We 
conclude  that  the  use  of  the  series  for  <p( c,  T)  can  greatly  facilitate  the  estimation  of  a  reliable 
value  for  T  from  a  calculation  of  the  first  few  terms  of  the  high-temperature  expansion  for  the 


K.  Dwight  H.E.  Stanley 

T.  A.  Kaplan  N.  Menyuk 
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TABLE  V-3 

RATIOS  OF  SUCCESSIVE  COEFFICIENTS 

IN  POWER  SERIES  REPRESENTATION 

OF  OPTIMIZED  CORRELATION  (S  =SB=3/2) 

u 

bA> 

b2Ai 

vs 

b4A3 

2.00 

1.000 

0.493 

0.606 

0.737 

1.60 

1.000 

0.627 

0.639 

0.623 

1.20 

1.000 

0.732 

0.690 

0.685 

0.80 

1.000 

0.810 

0.744 

0.754 

0.40 

1.000 

0.864 

0.787 

0.812 

0.00 

1.000 

0.899 

0.822 

0.835 

-0.40 

1.000 

0.917 

0.844 

0.850 

-0.80 

1.000 

0.925 

0.854 

0.858 

-1.20 

1.000 

0.924 

0.857 

0.856 

-1.60 

1.000 

0.918 

0.852 

0.849 

-2.00 

1.000 

0.910 

0.846 

0.839 

2.  Conduction-Electron  Spin  Density  Induced  by  Localized  Electrons 

4  5 

Recently,  Overhauser  and  Stearns  have  interpreted  Mossbauer  measurements  on  Fe-Al 

alloys.  They  concluded  that  there  is  a  large  qualitative  difference  between  the  true  4s -electron 
susceptibility  and  the  free-electron  susceptibility  Xq(<1K  due  to  interactions  between  4s  electrons. 
[The  free  electron  theory  gives  a  4s  spin  density  S(R)  which  is  in  reasonable  agreement  with 
experiment  for  lattice  points  R  ^  0,  but  which  is  about  an  order  of  magnitude  off  at  R  =  0.)  We 
have  reconsidered  the  theory  of  the  polarization  of  conduction  electrons(s)  by  exchange  inter¬ 
action  with  localized  electrons,  neglecting  s-s  interactions.  We  note  that  the  free-electron  ap¬ 
proximation  gives  conduction-electron  wave -functions  ip k  proportional  to  expfik  •  r]  at  the  very 
center  of  the  atomic  core  —  obviously  a  poor  approximation  for  interpreting  hyperfine  field 
measurements.  We  therefore  took  ip^  as  a  plane  wave  orthogonalized  to  the  core  for  determining 
the  value  of  ip ^  at  the  nuclear  positions.  We  found  that  although  S(R)  for  R  0  is  essentially 
unchanged  by  this  modification,  S(0)  is  changed  appreciably  to  be  more  in  accord  with  experi¬ 
ment.  These  calculations  show  clearly  that  the  conclusion  of  Overhauser  and  Stearns  was  pre¬ 
mature,  and  that  a  theory  neglecting  s-s  interactions  might  very  well  give  reasonably  good 
agreement  with  these  experiments.  ^  ^  Kaplan 


3.  Chemical  Inhomogeneities  and  Square  B-H  Loops 

A  paper  of  this  title  has  been  submitted  to  the  Journal  of  Applied  Physics.  The  abstract 
reads  as  follows: 
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"It  is  shown  that  in  ferrites  exhibiting  spontaneously  square  B-H  loops, 
reverse -domain  nucleation  may  occur  either  at  a  grain  boundary  or  at  a 
chemical  inhomogeneity  having  a  volume  roughly  0.1a  grain  volume,  a  mag¬ 
netization  about  2  percent  different  from  that  of  the  parent  matrix,  and  no 
static  crystallographic  discontinuity  at  the  inhomogeneity-matrix  interface. 

It  is  pointed  out  that  such  chemical  inhomogeneities  may  occur  in  spinels 
having  Jahn-Teller  ions  in  concentrations  that  are  too  small  to  produce  a 
cooperative,  static  distortion  of  the  entire  crystal  from  cubic  symmetry. 

It  is  postulated  that  the  formation  of  such  inhomogeneities  is  responsible 

for  the  sharp  knee  in  spontaneously  square  B-H  loops  found  in  many  ferrites 
3+  2+ 

containing  Mn  and  Cu  ions,  in  particular  the  magnesium-manganese 
composition  used  in  memory-core  ferrites." 

J.  B.  Good enough 

4.  Energy-Level  Scheme  for  CoTiO^  and  FeTiO^ 

2+  2  + 

Previous  DC  susceptibility  measurements  of  the  effective  moments  of  Co  and  Fe  ions 
in  CoTiO-  and  FeTiCL,  gave  values  of  5.47  and  5.21  Bohr  magnetons  for  the  respective  ions. 

Based  on  a  spin  of  3/2  for  Co  and  2  for  Fe  ,  the  spectroscopic  splitting  factors  necessary 
to  explain  these  moments  are  2.83  and  2.13,  respectively.  Since  the  effective  moment  of  CoTiO^ 
is  much  greater  than  the  spin-only  value  of  3.87  Bohr  magnetons,  there  should  be  a  substantial 
contribution  of  the  orbital  angular  momentum  to  the  moment,  while  in  FeTiO^  the  measured 
moment  exceeds  the  spin-only  value  (4.9  Bohr  magnetons)  by  a  much  smaller  amount  and  should 
have  a  smaller  orbital  angular  momentum  contribution. 

It  might  have  been  anticipated  from  the  anomalously  large  g-factor  in  CoTiO^  that  the  crystal- 
field  anisotropy  would  be  greater  in  CoTiO^  than  in  FeTiO^.  However,  this  is  inconsistent  with 
experiment.  The  dipole -dipole  anisotropy  favors  spin  alignment  within  the  basal,  c  planes.  In 
the  case  of  FeTiO^,  the  spins  are  directed  along  the  c-axis,  indicating  a  crystalline  anisotropy 
greater  than  the  dipole-dipole  anisotropy.  For  CoTiO^,  the  spins  lie  in  the  c-plane,  indicating 
a  crystalline-field  anisotropy  that  is  smaller  than  the  dipole -dipole  anisotropy. 

These  inconsistencies  are  partly  resolved  by  considering  the  ground  states  of  the  ions, 
which  in  an  octahedral  field  are  triply  degenerate,  as  being  split  under  the  action  of  spin-orbit 
coupling  into  a  doublet  ground  state  for  CoTiO^  (J*  =  l/2)  with  a  g*  =  (13/3);  and  a  triplet 
ground  state  for  FeTiO^  (J*  =  1),  with  a  g*  =  7/2  (Ref.  6).  In  CoTiO^,  this  will  lead  to  an  ef¬ 
fective  moment  which  will  be  particularly  temperature  dependent,  since  the  excited  states  having 

J*  =3/2  and  5/2  are  separated  from  the  ground  state  by  (3/4)  A  and  (2)  A,  respectively,  where 
2+  -1 

tie  Co  free-ion  value  of  A  is  equal  to  180  cm  .  A  calculation  of  the  effective  moments  at 
temperatures  of  60°  and  120  °K  gives  values  of  5.25  and  5.05  for  the  respective  moments,  show¬ 
ing  that  for  the  temperature  range  over  which  the  susceptibility  was  measured,  the  experimental 
moment  exceeds  the  computed  value  by  10  percent  at  most. 

The  measured  antiferromagnetic  resonance  frequency  f  for  CoTiO^,  extrapolated  to  0°K 
is  145kMcps.  With  octahedral  coordination  and  crystalline -field  contribution  to  the  anisotropy 
f:.eld  neglected  because  of  the  doublet  ground  state,  the  resonance  frequency  at  0°K  with  no  ap¬ 
plied  magnetic  field  is  given  by  f  =  e/2  me  (gV2>  —  K|),  where  K||  and  K|  are  the 
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Fig.  V-l .  Modified  ampoule  and  furnace 
w  temperature  profile  for  vapor  growth  of 
x  ReOg  crystals. 

z 


Fig.  V-2.  Temperature  dependence  of  resistivity 
for  an  ReC>3  single  crystal.  Measurements  made 
with  two  instruments  using  DC  currents  of  lOOma 
(•)  and  2  amp  (■). 
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anisotropy  energy  constants  for  spins  along  and  perpendicular  to  the  c-axis.  To  satisfy  the 
resonance  frequency  equation,  g*  must  equal  3.56,  which  differs  somewhat  from  the  octahedral 
field  value  of  13/3.  This  indicates  the  presence  of  a  trigonal -field  component  in  addition  to  the 
octahedral. 

A  trigonal  field  acting  on  Co  will  introduce  no  additional  anisotropy  in  the  ordered  state 
since  the  ground  state  which  is  well  separated  from  excited  states  has  J*  =  l/2.  However,  in 
FeTiO^,  where  the  ground  state  has  J*  =  1,  a  small  trigonal  field  lifts  the  degeneracy  only 
slightly  and  a  considerable  crystal  field  anisotropy  is  possible. 

J.  J.  Stickler 
J.  B.  Goodenough 
H.  J.  Zeiger 


B.  EXPERIMENTAL 

1.  Crystal  Growth  and  Electrical  Properties  of  Rhenium  Trioxide 

Interest  in  the  electrical  properties  of  rhenium  trioxide  has  been  indicated  in  previous  Solid 
7  8 

State  Research  Reports.  '  It  was  shown  that  this  material,  when  prepared  by  conventional  tech¬ 
niques,  is  usually  nonstoichiometric.  However,  effective  purification  was  obtained  via  a  vapor 
transport  technique  using  iodine  as  transporter.  In  these  experiments,  transport  of  ReO^  oc¬ 
curred  in  a  two-zone  furnace  under  conditions  of  high  iodine  vapor  pressure  in  the  direction  of 
the  cooler  zone.  These  observations  suggested  the  application  of  this  technique  to  the  growth 

of  ReO~  single  crystals.  Several  preliminary  experiments  were  based  on  the  method  described 
*  9 

by  Fischler  for  the  growth  of  tellurium-doped  silicon  crystals.  This  consists  of  growth  from 
the  vapor  by  slowly  pulling  a  sealed  ampoule  containing  the  charge  through  a  vertical  gradient. 
Although  these  experiments  indicated  that  the  technique  was  in  principle  applicable,  transport 
was  not  complete  and  the  crystals  obtained  were  too  small  for  conductivity  measurements.  These 
results  indicated  that  effective  transport  was  very  sensitive  to  the  temperature  differential  be¬ 
tween  the  reaction  and  crystallization  ends  of  the  ampoule  and  to  the  iodine  pressure.  Thus,  the 
rate  of  pull  was  too  rapid. 

In  order  to  establish  the  optimum  conditions  for  crystal  growth,  the  ampoule  has  been  modi¬ 
fied  as  shown  in  Fig.  V-l.  Platinum  and  platinum-10%  rhodium  thermocouple  leads  were  sealed 
into  the  top  of  the  ampoule  in  such  a  way  that  the  unconnected  ends  extended  about  1  mm  into  the 
ampoule  and  were  about  0.5  mm  apart.  Externally,  the  leads  were  connected  in  series  with  a 
flashlight  battery  and  bulb.  With  this  arrangement  the  ampoule  was  again  pulled  slowly  through 
the  temperature  gradient.  As  crystal  growth  began  at  the  top  of  the  ampoule,  the  deposited 
R.eO-j  shorted  the  internal  leads,  completing  the  circuit  as  indicated  by  the  glowing  bulb.  At  this 
point,  the  pulling  mechanism  was  stopped  and  the  growth  temperature  determined  by  connecting 
the  thermocouple  leads  to  a  potentiometer.  The  ampoule  was  maintained  in  this  position  for  a 
period  of  24  hours.  Under  these  conditions,  transport  was  complete  and  relatively  large  crystals 
(approximately  2.5  mm  on  an  edge)  suitable  for  conductivity  measurements  were  obtained.  The 
optimum  crystallization  temperature  determined  in  this  way  appeared  to  be  370°  ±  5  °C  with  a 
temperature  differential  at  the  bottom  of  the  ampoule  of  about  10°  to  15  °C. 

Electrical  resistivity  values  obtained  on  a  typical  ReO^  single  crystal  are  shown  in  Fig.  V-2 
as  a  function  of  temperature  between  77°  and  300°K.  These  values  are  remarkably  low  for 


59 


Section  V 


an  oxide.  The  room-temperature  resistivity  is  about  one  order  of  magnitude  lower  than  that  of 
the  most  highly  conducting  tungsten  bronze  and  the  77°  value  is  similar  to  that  of  electrolytically 
pure  silver.  In  accordance  with  the  premise  previously  stated/  these  results  for  ReO^  indicate 
that  occupancy  of  the  A-sites  in  oxides  with  perovskite  or  ReO^  structures  is  not  a  necessary 
prerequisite  for  the  formation  of  conduction  bands  in  these  materials. 

A.  Ferretti  R.W.  Germann 
D.B.  Rogers  E.J.  Delaney 

2.  The  Effect  of  Trivalent  Manganese  on  Crystal  Chemistry  of  Some  Lithium  Spinels 

The  following  is  the  abstract  of  an  article  that  has  been  submitted  to  the  Journal  of  Applied 
Physics: 

" Crystallographic  properties  of  the  spinel  systems  LiQ  ^Ga ^  ^ 

(0  <  x  <  0.7),  Li^  ^Fe^  5-x^nx^4^°  ^  x  <  0-5),  and  Lig  ^Al^ 

(0  <  x  <  0.5)  have  been  investigated  in  order  to  determine  the  effects  of  tri¬ 
valent  manganese  substitutions  on  the  crystal  chemistry  of  the  ordered  spinel 

+3 

hosts.  In  each  system  relatively  small  concentrations  of  Mn  removed  long- 
range  ordering  of  lithium  ions  on  the  octahedral  sites.  This  loss  of  order  was 
accompanied  by  anomalies  in  the  dependence  of  lattice  parameters  on  composi¬ 
tion.  These  observations  can  be  qualitatively  interpreted  by  assuming  that 
+3 

Mn  ions  tend  to  cluster  in  order  to  reduce  the  elastic  energies  associated 
with  Jahn-Teller  stabilizations.  These  effects  may  provide  indirect  evidence 
for  lattice  imperfections  such  as  have  been  invoked  to  account  for  square  B-H 
hysteresis  loops  in  ferrite  materials  containing  Jahn-Teller  ions.'1 

4-3 

In  the  course  of  this  discussion,  it  is  pointed  out  that  clusters  of  Mn  ions  such  as  those 

envisaged  are  consistent  with  the  mechanism  for  reverse  domain  nucleation  assumed  by  Good- 

enough  in  a  companion  paper.  D  B  Rogers 

R.  W.  Germann 
R.  J.  Arnott 


3.  Phase  Transitions  in  LaCoO^ 

LaCoO^  has  been  investigated  by  means  of  x-ray  powder  diffraction  in  the  temperature 
range  500 °K  <  T  <  1600 °K  and  by  means  of  differential  thermal  analysis  from  room  temperature 
to  1500°K.  The  two  methods  have  consistently  shown  the  existence  of  a  first-order  transition 
at  1210°K.  Electrical  measurements  performed  by  previous  authors  indicate  a  sharp  minimum 
in  the  resistivity  at  this  temperature.  Further  measurements  are  required  before  a  definitive 
interpretation  of  the  peculiar  magnetic,  electrical,  and  crystallographic  properties  of  this 
compound  can  be  established.  p  ^  Raccah 

C.  Anderson 

4.  Properties  of  B31  Compounds 

In  order  to  further  examine  the  dependence  of  the  signs  of  the  magnetic  couplings  and  the 
magnitude  of  the  atomic  moments  of  B31  compounds  as  a  function  of  the  electron/atom  ratio, 
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crystallographic  and  magnetic  studies  of  the  systems  Fe^Mn^^P  and  FeP-FeAs-FeSb  have  been 
initiated. 

Compounds  of  the  system  FexMn^_xP  were  prepared  by  heating  the  elements  in  a  sealed 
quartz  tube  to  1150°C  at  the  rate  of  20°<j/hour.  These  compounds  were  held  at  this  temperature 
for  three  days,  slow  cooled  to  800 °C  and  then  quenched.  Heating  below  800 °C  for  any  time  pro¬ 
duces  a  phase  separation.  The  unit  cell  decreases  slightly  with  increasing  x,  each  lattice  param¬ 
eter  obeying  Vegard's  law.  In  the  range  of  iron  Concentrations  0  <  x  <  0.25,  the  metamagnetic 
t  ransition  temperature  at  50 °K  in  MnP  increases  with  increasing  x  and  the  ferromagnetic  Curie 
temperature  T^  decreases  slightly.  The  paramagnetic  Curie  temperature  also  decreases  with 
increasing  x.  For  larger  x,  there  is  only  metamagnetism  below  the  long-range  magnetic- 
ordering  temperature,  the  paramagnetic  Curie  temperature  becoming  9  <  T^  and  continuing  to 
decrease  with  increasing  x.  The  Curie  temperature  becomes  negative  at  x  =  0.55  and  reaches 
a  minimum  at  x  «  0.8.  Iron  phosphate  is  ferromagnetic  with  ©  £  T  ,  suggestive  of  a  possible 
metamagnetic-ferromagnetic  two-phase  region  occurring  in  the  interval  0.8  <  x  <  1.0. 

J.  B.  Goodenough 
W.  A.  Newman 


5.  Temperature  Dependence  of  Attenuation  of  70-Gcps  Acoustic  Waves  in  Quartz 

With  present  techniques,  experimental  measurements  using  phonon  generation  at  70Gcps 
are  marginal.  Two  main  areas  are  now  being  investigated  that  may  lead  to  improved  system 
performance.  The  area  which  offers  the  greatest  possible  gain  is  improvement  of  the  over-all 
electromagnetic  to  acoustic  transducing  efficiency.  Thin-film  transducers  of  CdS  or  ferro¬ 
magnetic  films  with  a  magnetic  biasing  field  are  presently  available  alternatives  to  single-surface 
generation  of  the  quartz  itself.  Even  if  this  approach  does  not  yield  improved  efficiencies,  it 
permits  nonpiezoelectric  materials  to  be  measured.  Another  area  of  possible  improvement  is 
the  reduction  of  the  receiver  noise  figure  from  its  present  value  of  18  db  by  using  a  different 
mixer  configuration  or  by  preceding  the  mixer  with  a  klystron  preamplifier. 

The  present  system  is  also  being  modified  to  permit  attenuation  measurements  to  be  made 

downtol.5«K.  J.B.  Thaxter 

C.  D.  Parker 

P.  E.  Tannenwald 
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